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Abstract 
, ^ \ 

l 

Under NASA Grant NSG 1X88, .computer programs were developed for deter- 
mining constant output feedback gains > for linear systems subject* to both process 
and measurement uncertainties. However, in order to evaluate the effective- 
ness of these programs additional study was required of their applicability 
to the design of controllers for representative aircraft models. 

To this effect a sixth order linear longitudinal model and a 17 th order 
linear longitudinal model with five bending modes were used for the design of 
reduced state controllers for gust alleviation. 

Results show that the developed non-gradient based Zangwill-Powell op- 
timization program could indeed be used to design satisfactory output feedback 
controllers taking into account the needs for reducing vertical accelerations 
and structural loading effects i 
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1. INTRODUCTION 

Because application of conventional optimal linear regulator theory to 
flight control systems requires the capability for measuring and feeding hack 
the entire state vector, it is of interest to consider procedures for comput- 
ing controls which are restricted to be linear feedback functions of a lower 
dimensional output vector. Such a procedure, however, has its limitations in 
that the feedback gains will be functions of the initial state vector. In 
addition, the presence of measurement noise and process uncertainty can lead 
to additional problems relating to both modelling and computation. 

In order that such reduced state feedback control laws might be consid- 
ered for the design of flight control systems, an extensive study effort was 
conducted between 1 June 1975 and 30 Nov. 1976 under NASA Grant NSG 1188 mon- 
itored by Mr. Ray Hood of the Langley Research Center 

The objectives of this grant were to: 

• Develop procedures for computing optimal constant out put feedback 
gains for linear flight control systems taking into account the 
presence of turbulence, pilot commands, parameter uncertainty, 
and measurement noise. 

• Implement, on NASA Langley’s CDC6600, computer programs capable 
of determining such optimal output feedback gains. 

• Compare the performance of the various computational algorithms 
and investigate various procedures for modelling the system 
uncertainties. 

a Document the program operation for public dissemination. 

Toward these goals, working programs were developed for finding output 
feedback gains for linear systems in order to minimize both infinite and finite 
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time performance indices; ' Computational procedures included algorithms 
which require direct computation of the gradient of the index with respect 

h) . . (M 

to the gains and the algorithm proposed by Powell and modified by Zangwill 
which does not require gradient computation. Two different procedures for 
modelling plant parameter uncertainties were considered. 

Significant among the conclusions resulting from these studies were 
the following: ^ 

• Use of a finite time performance index can result in a set of gains 
which do not stabilize the closed loop system. 

e If it is possible at all to stabilize the system with the specified 
feedback configuration, then the optimization of an infinite time 
performance index will yield a set of gains that do indeed stabilize 
the closed loop system. 

• Optimization of infinite time indices is less time consuming than 
the optimization of finite time indices because of the need to solve 
algebraic and not differential equations. However, the computation 
procedures for infinite time indices do require initialization with 
a gain matrix that stabilizes the closed loop system. 

« Application of 'gradient based algorithms to the optimization of 

infinite time indices can result in divergence. This in particular 
results when the gradient is computed using the algebraic steady 
state Riccati solutions to the matrix covariance and co-state equa- 
tions. These solutions are only valid provided that the gain matrix 
under consideration stabilizes the closed loop system. Consequently, 
if during the search process a perturbation results in a destabili- 
zing gain, then the corresponding computed gradient will be meaningless. 
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Although these previous efforts resulted in a set of computer programs 
which can be used for finding a set of gains for a given reduced state feed- 
back control configuration, it was necessary that further explorations be made 
of their utility to flight control system design. In particular the effects 
of modelling, sensitivity, and stabilization needed consideration with respect 
to more realistic aircraft models. 

Thus under the present grant, IfSG 138^, further studies have been made 
towards the application of the reduced state stochastic infinite time optimi- 
zation programs to the design of control systems for representative flexible 
aircraft. Towards this goal the following tasks have been performed: 

• Incorporation of a procedure which computes if possible an initial 
gain matrix. 

• Further comparison of both gradient and non-gradient based pro- 
cedures for designing reduced state feedback flight control systems. 

• Evaluation of the reduced state feedback control computation pack- 
age in designing a gust alleviation controller for a representative 
flexible aircraft. 

Results show that the non-gradient based reduced state feedback control design 
program can indeed be used for designing acceptable controllers for a 17th 
order flexible aircraft model. This conclusion was based upon the perfor- 
mance of controllers designed for the reduction of vertical acceleration and 
structural loads in the presence of a vertical wind gust. 
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2. Problem Statement 
2,1 Process Model 

The optimization algorithms discussed in ref. 1 are applicable to systems 
described by the following set of state equations: 


Process : 
Reference: 
Disturbance: 
Where: 


X = 


X = 


X = 
n 


X = 


X = 


X = 
n 

U = 


= A X 

+ B 

U + G 

P P 

P 

P P 

= A X 

+ B 

w 

r r 

r 

r 

= A X 

+ B 

w 

n n 

n 

n 

= (WXP 

x l) 

plant 


n 


(WXR'x l) reference state 
(WXET x 1.) disturbance state 
(WUP x 1) control vector 


( 1 ) 

( 2 ) 

( 3 ) 


Wp = (NXP x l) . process white noise vector with covariance 


w 


= (WWR x l) reference white noise input with covariance W, 


R 


w n = (WWW x 1) disturbance white noise input with covariance W q 


Given the system defined by equations 1, 2, and 3> the available outputs 
are to be designated as: 


Y = C X+C X+Y 
p pp p pn n 


Y = C X 
r rr r 


Y - C X 
n nn n 


w 

( 5 ) 

(6) 


where Y^ = (WYP x l) process measurement 

Y = (WYP x l) reference measurement 


Y^ = (WYW x 1) disturbance measurement 


and y = a zero mean white noise disturbance with covariance I\ 
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The control U is to be of the form: 
P 


U=K Y + K Y + K, Y 
» y D p y r r y n n 


(7) 


where the gain matrices K . K ? K are to be computed so as to minimize: 

y y y 

J r J n 


J = Lim £,[ ~ f {(Z - Z ) T §(Z - Z ) + TJ T R U } dt] 
j j n « r' p r' p p’ 


( 8 ) 




and where £, [• ] denotes statistical expectation and the controlled variables 

: Z = D X 


Z and Z are of the form: Z = D X + D X 

p r P PP P pn n 


rr r 


2.2 Optimization Procedures 

For notational convenience, eqs. 1, 2, 3 will be compressed into the form: 
X = A X + B U + V ± (t) (9) 


I A 


where X - 


X 

E 

X 


VI 

\(t) = /w p + V(t)\ 


B W 
r r 


i B W i 
\ n "■ / 


£ (v^T) = V 
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p 
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B = 

(M 

0 

A 

r 

0 


0 

i 


0 



(° i 


Similarly, eqs. k, 5, 6 can be expressed as: 

Y = C X + n 
and the controlled variables as: 

Z = DX 


(10a) 


(10b) 
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where Y = 


M 

n = 

Y \ 

c = 

h 

0 

C \ 

D «/ 

p \ 


\ 


PP 


pn \ 

y r 


° 


0 

c 

0 


U, 


' 0 / 


\ 0 

rr 

0 

C 1 


\ 1 



\ 


nn / \ 


D 'O D 

PP 

0 D 0 

rr 


pn 


z - (z p , Z r ) 

Therefore the control TJ can be rewritten as: 

P 


TJ = KY 

P 


(11) 


where K = [K , K , K ] 


or 


V V y n 


With this notation, eq. 8 for the performance index can he rewritten as: 

T 

J = Lim £ [ ~ f [Z T q Z + u T r u ] dt] 

3. a p p 


T-*oo 


where Q = 
a 


Q -Q, 
-Q 
0 


Q 0 
0 0 


T-*x> 


J = Lim £ J [}? D T Q a D X + X 1 C T H? R K 'C 

+ 2 5? C T K 1 E n + n T K T R K n] dt] 


•c X 


Note that some simplification is possible since 
£ (xWR n) = 0 


Thus, the actual index selected for minimization was: 
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T 

J = IAm £. ~ f X T (D T Q D + C T K T R K C) X dt 
_ 1 ^ ry a 

I X L ^oo 0 

+ £ (n T K T R K n) 


or^ 8,9 ^ J = Trace [(D T U + C T K T R K C) P ] + Trace (K T R K R) (12) 


where R = covariance of n 


T 


and P g = Lim “ J £ (X X T ) 


dt 


Q>*co 


Rote 1 that the last term in eq. 12 is a penalty on the gains. Increasing R 
and/or R therefore has the effect of reducing the magnitude of the resulting 
gains. If this integral does not converge, then it is common to use the 
steady state value . 


P a = Lim (X(t) X T (t)). 


It is of importance to note that provided (a_b k C) is a stable matrix, 


may 


/ON 

be found by solving the linear matrix equation: ' ' 


0 = (A-B K C) P + P (A-B K C) T + B K R K T B T + V (13) 

- S S 


Thus for a given set of gains it is possible to compute the performance 

index from ( 13 ) if the steady state covariance matrix P_ is available. 

b 
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3. Computational Procedures 

3*1 Non-gradient based Optimization 

Since the performance index (eq. 12) is easily evaluated given a value 

(o 4) 

for the gain matrix K, the Zangwill-Powell’ ' method which does not require 
gradient computation is very attractive for optimization. In particular, the 
IMSL sub-routine ZXPCWL was used for implementation. 

Letting UK denote the number of unknown gains to be determined, each 
iteration of the procedure begins with NK single dimension minimizations 
along NK linearly independent directions . . . , each starting from 

K q , the most recently available gain matrix. Initially K q is user selected 
and the directions correspond to the coordinate vectors so that during the 


first set of minimizations only one gain element is changed at a time. Sub- 
sequent iterations define the direction vector set (§g, ..., 5^, §), 

where § is chosen such if the objective were quadratic, after k iterations, 


the last k of the direction vectors would be mutually conjugate. These re- 

(3 4) 

vised directions are then used for the next iteration. ’ ' 


This procedure is especially useful for the index defined by eqs. 12, 13- 
In this case if the algebraic solution of the steady state version of eq. 13 
is to be meaningful, then the gain K must stabilize the closed loop system 
(A-B K C). If K is not stabilizing, then J = o. If during the process .of 


searching along a particular direction vector 5 ^, perturbations are such that 
(A-B K C) becomes unstable, then under program control, the size of the per- 
turbation will be reduced (to zero if necessary) prior to the subsequent search 


along the next direction vector. This is in contrast with conventional grad- 


ient type search procedures which do not have other search direction vectors 

(5) 

available when instability results. 
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To account -for the above stability problem using the program ZXPOWL, 

the eigenvalues of (A-B K C) were computed for each perturbed value of K, 

50 

and J was then set equal to 10 whenever an eigenvalue was found to be greater 
than or equal to zero. 

To insure convergence of the procedure, it is necessary that the initial 
gain K° stabilize the closed loop matrix, (A-B K°C). Such a matrix, if 
one exists at all for the permissible feedback structure, may be found by 
defining an initial phase to the procedure in which the performance index 
to be minimized with respect to K° is 

= Maximum real part of the eigenvalues of (A-B K°C). 

As soon as a gain vector is found which' is such that all eigenvalues of (A-B K°C) 
have negative real parts, then the final phase can be undertaken in which the 
index is that defined by eqs. 12 , 13 . 

3-2 Gradient based computation procedures 

In order to compare the performance of the non-gradient based Zangwill- 
Powell procedures with gradient based algorithms, a program implementing the 
sequential unconstrained minimization technique (SUMT)^^ was also considered. 

This program is modular in structure so as to facilitate changes, and 
uses a series of control, special purpose, and user supplied subroutines in 
order to solve the general mathematical programming problem: 

Determine the vector X so as to 
(a) minimize F(X) 
subject to: 

R.(X) >0 j = 1 , ..., M 

t) 

R.(X) = 0 
3 - 


3 = M + 1, .. ., (M + MZ) 
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This is done by solving a sequence of unconstrained problems whose 
solutions approach the solution of (a). 

Specifically , the SUMT procedure uses the function defined as: 


M+MZ 


AM ■] AM « A’JiJ 

P(X,r) = F(X) - r S ln(R (X)) + ± Z (E (X)) 


j=M+l 


Using a designated search procedure, a sequence of, P(X,r.) is minimized 


for r = x v r. 


r, , where r . , = r. / and c > 1. Under suitable condi - 

K 1+1 1/C 


tions, the minima of P represented by X( r ]_)j •••» X(r^) exist and 

approach a solution of the constrained problem (a) as r fe -» 0. ^ That is: 

Lim X( r v )) = X* 

r: -0 K 
b 


Lim F(x(r )) = F(X ) 

r, -0 
k 


Note that the equality constraints are satisfied only in the limit as 
r^"*0 - An extrapolation option is available which can, in some cases, accel- 
erate the convergence. To start with, initial values X for X and r, for r, 

—o—l 

must be available. These can either both be read into the program, or r^ 
can be computed using one of two options which have been programmed; namely: 


OPTION 1: Find r to minimize 


P <^> rh' 1 \ p (^> *> 


This is useful if at least one 


R.(X ) « 0 
2 ~o 


OPTION 2: Ignore the equality constraints and minimize 


| |V P(X , r) 
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Furthermore, if results in one or more of the inequality constraints not 
being satisfied, then the program operates in a feasibility phase by forming 
an auxiliary objective function equal to the negative of the sum of all the 
violated inequality constraints. When a constraint is noted to become feasible 
during the minimization of the auxiliary objective, it is removed and included 
in the effective constraint set. 

To apply the StfMT program to the design of a stochastic reduced state 

feedback controller, the vector of unknowns x must be associated with the 

elements in the gain matrix and the constraints R. must be formulated so as 

3 

to yield a stable closed loop system. In particular, the equality constraints 
were used, and R^ was formulated to be the negative of the maximum real part 
over all eigenvalues of the closed loop system matrix. 
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4. Results and Discussions 

Although most of the contractual effort was expended on the application of 
optimal reduced state feedback gains to realistic aircraft models , some prelim- 
inary activities were devoted to finding procedures for computing an initial 
stabilizing gain matrix and to comparing the gradient based SUMT procedures: with 
the Zangwill-Powell approach. 

4.1 Gain initialization 

Since lose of the Zangwill-Powell method in conjunction with eqs. 12,13 

requires that the initial gain matrix stabilize the closed loop matrix (A-B K C), 

it is important that a procedure be incorporated for finding such a gain matrix, 

if it exists, for a specified feedback structure. 

(12) 

After studying several possibilities, it was decided to use the Zangwill- 
Powell procedure to minimize , with respect the gains , the maximum real part of 
the closed loop eigenvalues. This procedure has been incorporated into the 
program, and is called, if the original specified gain yields an eigenvalue with 
a positive real part. If after a specified number of iterations (ITMAX), the 
Zangwill-Powell procedure does not find a stabilizing gain matrix, then a message 
is printed, and the program stops. 

4.2 Gradient based optimization procedures 

In order to compare the performance of the non-gradient based Zangwill- 

Powell procedures with gradient based algorithms, a program implementing the 

sequential unconstrained minimization technique (SUMT) discussed in Section 3-2 

was developed and applied to both third and sixth order linearized longitudinal 
(13) 

models. In general, it was observed that the SUMT procedure required more 

computer time than the Zangwill- Powell method to converge; m fact, in many cases 

(13) 

the final SUMT performance index exceeded that reached by the Z-P procedure. 

These relative inefficiencies were attributed not only to the need for gradient 
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computation, but also to the requirement for optimizing a sequence of many 
unconstrained problems. 

4.3 Application to representative aircraft models 

In order to evaluate the effectiveness of the infinite time reduced 
state feedback controller program, various controller configurations were 
developed according to NASA suggestions on two analytical models of the TIFS 
aircraft. The TIFS is a fly by wire C-131 aircraft owned by CALSPAN. 

With its onboard computer and separate controller for all six rigid body 
degrees of freedom, it is a unique facility for control system and handling 
qualities research. 


4.3.1 Gust alleviation using a sixth order longitudinal model 
4. 3- 1.1 Flight control problem definition 

Initially a modified six-dimensional version of the TIFS aircraft ^per- 
turbed by a vertical wind gust was used for evaluation. The correspond- 
ing variable definitions in accordance with eqs. 1-6 were as follows: 
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Controlled outputs: 


( n zl\ _ j vertical acceleration at point $ 1 \ 

z ^l l vertical acceleration at point # 2 J 

The structural matrices (eqs. 1-6) used for design corresponded to the climb 
condition, i.e., h = 1524 m, V = ICo m/s, and are (units in radians and ft/sec.) 


A “ 

TO 




. 9 

A 0 

AV 

Aa 

6 e ' 

5 c 

—.1686 

.000035 

.000231 

VO 

CO 

-3* 

* 

1 

-4.3773 

-.19948 

1. 

0 . 

0 . 

0. 

0. 

0 . 

0 . 

-32.17 

-.0143 

18.027 

0. 

-3.0933 

1 . 

0.000013 

-.000531 

-1.223 

-.1273 

-.2667 

0. 

0. 

0. 

0. 

-20. 

0 . 

0 . 

0.- 

0 . 

0. 

0. 

-4o. 

0 . 

0. 



.486 

t 



G - 
P 


0 

18.027 

-1.223 

0. 

0. 


An = -.2784 , Bn = .01815 


°PP 


.( 64.63 

■^pp — \ -61.82 

_ /444.2\ 
P n “ \407.8 j 


.00318 

.00580 


.176 

.193 


444.2 

407.8 


212.1 

-116.2 


100.4 \ 
85.5 J 


Sensor noise standard deviations: 

a - .5 deg/sec 
9 

0 Q = .2 deg. 

0 = 0 - .2 deg. 

a ag 



For design purposes the gust was selected to correspond with a standard 
deviation vertical wind speed of 15 feet/second. However because of a pro- 
gramming error (discussed in more detail in 4.3.2) the computed gust standard 
deviation was in reality 15/*/*, or 8.46 ft/sec = 2.58 m/s. Thus the gust 
would correspond to a situation encountered somewhere between normal and cumulus 
type sky conditions. 

The performance index weights were then chosen to be approximately equal 
to the inverse square of the maximum permissible values of the associated 
variables. This resulted in 


T 


J = lim 
T-® 


1 

T 


E 


J Hg 

J o 16 2 


n 32 2 

16 2 


6 2 
ec 

.4 2 


6 2 
ZC V 

,6 2 ’ 


dt 


(1U) 


4. 3- 1. 2 Simulation Results 

Initial studies concerned with the effects of various feedback signals 
were followed up with a study' of the influence of sensor noise and different 
initial gains. Then using a typical feedback configuration, the sensitivity 
to flight condition changes and gust variance was also investigated. In each 
case evaluations were made relative to the resulting value of the performance 
index and the reduction in root mean square vertical acceleration, measured 
over five seconds; comparisons were made with respect to the open loop values: 

n zl = 9-75 ft/sec. 2 

n^ - 9*28 ft/sec. 2 

Effect of feedback configuration 

The following combination of sensed outputs were e xami ned with regard 
to their effectiveness in gust alleviation: 

* 

The same pseudo white noise sequence w q was used in all evaluations. 
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q, 0, a 
q, 0, (a-ag) 
q, 0, a, ag 
q, 6, (owxg), ag 

Examination of the results shown in Table 1 corresponding to an open loop or 
zero gain initialization show that all configurations reduced the rms vertical 
acceleration better than 60 $ and that use of the gust induced angle of attack 
ag in the feedback (either directly or implicitly through a-ag) further im- 
proved this effectiveness by an additional 20-40$. Note also that in the 
presence of sensor noise with the intensities cited in Section 4.3.1, the rms 
values of the vertical acceleration became moderately large. This motivated 
an accounting of the sensor noise characteristics in the optimization formula- 
tion. 

Effects of sensor noise 

Based upon the above results, the sensor noise was accounted for by 
assigning the appropriate variances to the noise term y of eq. 4. Related 
results shown in Table 2 show that taking into account the, sensor noise in 
the design phase results in the presence of distinct short period and phugoid 
modes along with smaller gains and a tremendous reduction (with respect to Table l) 
in the vertical acceleration responses. Note, however, that in order to 
achieve any significant improvement over the open loop case, it becomes neces- 
sary to include either a direct or implicit feedback of the gust effect ag. 

Sensitivity Studies 

In order to evaluate the sensitivity of the design procedures to varia- 
tions in the initialization and in the model, various perturbations were made 
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to the scenarios of the preceding two cases. First to assess the effect of 
different initial gains, a gain K° was found (as discussed in section 4. 1 to 
minimize the maximum real part of the eigenvalues of (A-BK°C). This gave 

Q / -.001616 -.09349 -.002438 \ 

\ -.0000851 -.002114 .0003967 J 

with eigenvalues: - 20.01 

-4o. 

- 1.345 ±3 .7702 

- .1090 ±3 .00412 

- .2784 

Using this gain to initialize the optimization in the presence of sensor 
noise led to results almost identical with those presented in Table 2. How- 
ever, without sensor noise, the optimal feedback gain between 6 and q was 

z 

in all cases greater than 1500 (in magnitude). 

To test the sensitivity of the gains with respect to flight condition 
changes, the Table 1 gains for the (q, 0, a-ag) feedback configuration which 
was derived for the climb condition was used in simulations at landing 
(h=6lm, ¥=68 m/s)and cruise (h-3048m, V=150m/s). This resulted in the follow- 
ing rms errors in vertical acceleration: 



n 31 

n 32 

Landing 

6.20 

4.92 

Cruise 

2.85 

2.80 


It was noted that although these values were significantly lower than the 

corresponding open loop values: 

n 31 n 32 

Landing 17.05 12.95 

Cruise 12.73 12.25- 

* they were higher than the corresponding Table 1 values for the climb conditions. 
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Thus, the gains for the (q, 0, a-ctg) configurations were redesigned for the 
climb conditions “but with a process noise term included to account for the 
uncertainty. To reflect the expected degree of uncertainty, the covariance 
of this fictitious process noise was chosen as: 

DIAG (.00006, 0., .00003, .00003, 0.,0.) 

This resulted in a set of smaller gains which gave the following rms vertical 
accelerations: 



n 31 

*32 

landing 

5.31 

3.27 

Climb 

2.20 

2.20 

Cruise 

3.02 

3-05 


Note that the inclusion of the process uncertainty improved the response at 
the landing configuration at the expense of small degradations in the response 
at the other two flight conditions. 

' Finally, in order to determine the effect of the assigned gust variance 
on the performance, gains were redesigned for the (a, 0, C£-Q!g) feedback: con- 
figuration with the assigned gust variance increased by a factor of four (rms 
value = 30 fb/sec). Although the gains were approximately doubled in value, 

it was determined through simulation -that both the larger gust gains and the 
smaller gust gains improved the open loop response by about a factor of three 
regardless of which gust was being applied. 

4.3.1. 3 Discussion 

- 

On the basis of the above results for the sixth order system, it was 
concluded that stochastic reduced state feedback design procedures would be 
useful in the design of aircraft gust alleviation control systems. Thus con- 
sideration was subsequently given to high order models which take into account 
bending modes, additional control surface deflections, and loading effects. 
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4.3-2 Reduced State Feedback Control of a flexible aircraft 
4. 3.2.1 Problem Definition 

To more realistically evaluate the use of the reduced state feedback 
control program, consideration was directed towards a 17th order model of the 
TIFS aircraft which incorporated five bending modes and three first order 

( 14 } 

actuators. ’ This model for the unaugmented TIFS was derived by CALSPAN 

using a quasistatic reduction on the equations which had been obtained with 
the FLEXSTAB estimation program. 

This data taken from Appendix B of ref. 14 (and reproduced in Appendix B of this 
report) for each of the 3 flight conditions (cruise, 'climb, land) was arranged into 

(15) 


the format: 


state X = A X +BU + G X 
eqn: p p p p p p n 


gust 

eqn: 


X = A X + B W 
n n n n n 


observation: 


Y = C X + C X 
p pp p pn n 


(16) 


(17) 


control: 


U = KY 
P P 


( 18 ) 


X^ is the aircraft state vector in the body axis system with components: 
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u 

x- velocity (m/s) 

w 

z-velocity (m/s) 

q. 

pitch r.ate (r/s) 

0 

pitch (r) 

n i 

1st bending mode 

“i 

i 

n 2 

2nd bending mode 

n 2 

n 3 

3rd 'bending mode 

"3 

n 4 

4th bending mode 

% 

n 5 

5th bending mode 

K 

6 

sa 

•symmetric aileron (deg) 

6 

z 

direct lift flap (deg) 

6 

e 

elevator (deg) 
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Note that the control surface deflections (all in deg) are treated as 

states from the actuator dynamics which are forced by the computed command 

U = (6 ,6,6). 

p sa 3 z e 
c c c 

The observation vector y consists of the components: 


zcg 

n zt 

n zwt 

n zFSF 

zast 

n zRHT' 



t r 


vertical acceleration (VA) at pilot station (g) 
VA at eg (g) 

VA at tail (g) 

VA at wing tip (g) 

VA on side force surface (g) 

VA at point on tail (g) 

VA at right horizontal tail (g) 
pitch rate at eg (deg/ s) 

Vane angle of attack (deg) _ 

Root shear (n) 

Root Bending Moment (n-m) 

/ 

i 

Root Torque (n-m) 
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The vertical wind gust x^, was generated according to the following 
equation received from NASA: 



where V = trim velocity 
o 

L = effective length 

and ^(t) is a zero mean, unit variance white noise process 


In actuality since the steady state variance of X resulting from this 

P n 

(7 2 
equation is selection of parameters so as to yield a specified c will 

really result in a variance lower by a factor of it* Consequently in the 

ensuing results which were derived for a " lm/s^ the true gust variance was 

o p 

actually .3183' m./s . Thus in order to extrapolate the normalized results to 
typical sky conditions, the following multiplication factors should bo used: 


Sky condition 

Normal ( 0 ^. « 6 ft/s) 
n 

Cumulus (a x ** 15 ft/s) 
n 

Thunderstorm (a^ » 3® ft/s) 
n 


Factor 

3 

8 

16 


From conversations with personnel at NASA Langley, it was determined 
that it would be of interest to determine reduced state feedback control gains 
for minimizing vertical acceleration and wing root bending moment which result 
from a vertical wind gust. To this effect the following two performance indices 
were considered: 

T 

J, = Lim ~ P f Flo 6 (BMR) 2 + 10" 6 [ 46 0 2 + 66 2 + 256 2 ]1 dt and (19) 

1 L J o L a c zc ec J 

T 

J = Lim | P f M T) 2 + 10 6 ri 2 ] + 10~ 6 [46 2 + 66 2 + 256 2 f| dt (20) 

2 T 0 J Q L zcg zp sa c ZC ec J 



Weights were selected so as to result in significant reduction (from 
open loop) in vertical acceleration and wing root "bending moment, without 
having the extrapolated control values for large gusts exceed the following 
limits: 


5 

< 

30° 

sa 

6 

< 

30° 

z 

5 

e 

< 

h- 1 
O 
O 


Results were evaluated on the "basis of a five second simulation run 
and upon the computed steady state covariances of the states and the penalized 
outputs . 


h. 3.2.2 Simulation results 

Evaluation of the reduced state feedback control design procedures 
consisted of a series of simulation studies which considered: 

• Initialization procedures 

• Comparison of various feedback configurations including full state 
feedback controllers computed using linear optimal quadratic regu- 
lator theory. 

• Effects of sensor noise 

• Treatment of the C response 

• Sensitivity to flight condition changes 
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Initialization Procedures 

Because several local minimum values for the gain matrix can result 

when using reduced state feedback gain optimisation procedures, it was of 

interest to examine various initialization procedures. To this effect the 

index ^(penalizing vertical accelerations) was minimized for the feedback 
T 

configuration y = (u, '0, n , q , 0! , to ). Initial gains were selected as: 

3 3 zcg ? -eg 3 v 7 g 

• All zero (i.e., open loop) 

• Those which resulted from optimization of the configurations 

T 

y = (n^ , l C gj a Y > m g) smarting from open loop. 

o Those which resulted from optimization of the configurations 
T 

y=(n , q ,a,oj,u, n.) 

J v zcg 3 ^cg’ v 3 g zwt ' 

with the gains on u and set to zero. 

Results shown in Table 3 indicate that it might, be desirable to optimize 
from more than one initial guess an'd to select the most satisfactory set of 
gains. Consequently it was concluded that a reasonable procedure for optimizing 
the gains for some specified configuration might be xhe following: 

1) First optimize the configuration consisting of only one output . 

2) Use the resulting gain from step one to initialize the optimization 
of a configuration consisting of two outputs. 

3) Continue adding one more output until the desired configuration 
has been optimized. 

k) If desired, re-optimize this configuration with all gains initialized 


at zero. 
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Effects of Feedback Configuration 

In order to examine the influence of a changing feedback structure, 

the two indices and defined by eqs. 19 * 20 were used for designing 

controller based upon feedback, of selected combinations ofn , n . . a , 

zcg 3 zwt’ -eg* 

- * 
a , cd , u, 0, and BMR. 

* o 

Controllers were designed based upon both an open loop initialization 
and upon the previously recommended procedure which involves the optimiza- 
tion of a sequence of an increasing number of outputs. In terms of convergence, 
performance index, and gust alleviation properties, the latter approach was 
superior. Corresponding results presented in Table 4 for designs based on 
minimization of J^, which in essence, penalizes vertical acceleration, indicate: 


• Gust feedback is desirable for reducing the acceleration. 

• Feedback of n .in addition to n , cd , q - and a significantly 

zwt zcg 5 g’ H cg’ v & ^ 

reduces both the vertical accelerations and the wing root bending- 
moment . 


© 


o 


Little is gained and at the expense- of increased controls by including the 

feedback of either 0, u, or BMR to the configuration defined by 

n , n . , cd , q , a . 
zcg’ zwt’ g 3 ^-cg 3 v 

Configurations including n^^ feedback show, compared with open loop, 
an increased damping of the lowest bending mode and an increase in 
the short period natural frequency. 


Results shown in Table 5 for configurations designed to minimize wing 
root bending moment indicate: 


p 

Note that in Table 2 and 3,™^ has been scaled by 10 , and BMR has been 

scaled by 10 . This was necessary since the initial gain perturbations 

performed by the program were too large in magnitude in both the positive 
and negative directions to indicate a meani gful search direction. 



• Feedback of n , cd , q , and a is almost equally effective as 
zcg’ g’ ^cg’ v- * ^ 

feedback of these quantities plus either u or 8 or BMR. 

9 Relative small decreases in BMR sometimes are accompanied by extremely 

high increases in n ‘and n 

zcg zp 

» The lowest order bending mode becomes more damped with feedback of 

n , id , q , and a ; however a distinct short period set of poles 
zcg g eg V 

is not distinguishable. 

.Simultaneous comparison of Tables h and 5 show that; 

9 Gains developed by penalizing BMR only (J^) can result in a BMR co- 
variance of about an order of magnitude less than that which results 
from penalizing the vertical accelerations (J^). However the corres- 
ponding shear and vertical acceleration covariances (from penalizing 
J^) are much larger than, those corresponding to open loop operation. 

9 Penalizing n and n (j_) results in vertical accelerations which 
are an order of magnitude or more smaller than those corresponding to 
open loop operation, accompanied by some reduction in both BMR and S„. 
9 Controls designed for penalizing BMR are larger than controls designed 
for penalizing the vertical acceleration. 

With regard to timing information, the number of iterations required 
by the Zangwill-Powell procedure are shown in Table ha. For 7 feedback out- 
puts (i.e., 21 gains) the computational time per iteration on the CDC-7600 
was about 2 minutes, for 3 feedbacks the time was about 1 -min per iteration 
and for one feedback the time was- about l/2 minute per iteration. Thus it 
was not unusual to use over an hour of computing time to design a 7 output 
controller . 

Full State Feedback Controller Comparison 
As a further evaluation of the effectiveness of reduced state feedback, 
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a comparison was made (using the 17 th order model) with controllers designed 
using full state feedback, i.e., u = KX. 

Note that included in this control would be direct feedback from the bending 
modes and actuater. 

Initially the same Q and R which define were used, but this resulted 
in extremely large controls; consequently R was increased until the control 
covariances were the same order of magnitude as those in Table 2. Note for 
example that the first entry of Table 6 for an R of DIAG (4,6,25) x 10 
results in extremely low vertical accelerations but at the expense of un- 
acceptable control magnitude. 

For the case R = DIAG (4,6,25), it should be noted that the correspond- 

t 

ing vertical acceleration, bending moment and control covariance are within 

an order of magnitude of those found in Table 4 for the feedback of n , cn , 

0 zcg g' 

Q 3 cc , n . . Thus a reduced state feedback controller if properly designed 
^cg 5 v* zwt 

can be almost as effective as a more complex and often unrealizable full state 
feedback configuration. 

With regard to control signal magnitudes, the gains developed from 

T 

minimizing J g for the configuration y = ( n zcg > <u g , n zw t’ S ave 

the following maximum values over a 5 second simulation for a gust of variance 
O .3183 m 2 /s 2 . 


max value in deg 


5 

sa 

6 

z 


These values except for perhaps 6 
storm situation (multiples « 16 ). 


I .38 

8.12 

.511 


should be acceptable even for a thunder- 

Further reduction in 6 would be achievable 

z 


by additional weighting in the performance index. 
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Effects of sensor noise 

Tables "J& and 7b depict the effects of modelling angle of attack and 
pitch rate sensor noise in the design phase for the output vector y - (n 

zcg 

T 

q , a , a) ) , In each case, the optimization procedure was initialized at 
eg v g 

open loop. 

To be noted are the following results: 


• Modelling of the sensor noise results in a considerable reduction 
in the gains which multiply q and a^. 

• As with Phe sixth order model, inclusion of sensor noise results in 

both a distinct short period and a distinct phugoid mode in the closed loop 
system. 

o Gains designed with sensor noise modelled, when compared with gains 
not incorporating sensor noise, gave a 3 - ^ order of magnitude reduc- 
tion in the steady state covariances of n , n , BMR, S . These 
a zcg 5 zp 5 R 

covariances were compared assuming a vertical wind gust input and 

the feedback of the a and q sensor noise components. 

v ^cg 


Design for C response 

For additional evaluation of flight control design using the reduced 

state feedback computation procedures, consideration was given to the C 

(l6) * 

response characteristic. ' Typically a step C would be commanded and 
both feedforward and feedback gains would be designed so as to yield accept- 
able transient and steady state specifications. However since the gain com- 
putation programs require asymptotic stability of the augmented process 
(including any reference system), step inputs could not be directly accomo- 
dated in the design. Instead, gains were determined so as to transfer the 
C response from some non-zero initial value to a near zero final value. 
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The resulting response transient characteristics are then identical to those 

that would result if the objective were to regulate the difference between 
* 

the actual C and some attainable steady state value. 

Using the relationship 

C* = 400q + n (4) 

eg zp 

where is in rad/sec and n^ is in g's, gains multiplying n ZC g> and 

a were determined so as to minimize 
v 

CO 

J, = if tA(C*) 2 + IQ- 6 [46 3a / + S6 sc 2 + 256 ac 2 ]] at 
o 

6 8 

for which Q = 10 and Q = 10 . Evaluation of the gains - was based primarily 

* 

upon a 5 second response to a unit C initial condiii on in the absence of 
measurement noise and disturbances. 

Results shown in Table 8 and in fig. 1 indicate: 

* 

e The C response resulting from the reduced state feedback design 

settles out in about half the time and with about half the overshoot 


which result under open loop control. 

g 8 

e Both weights (Q = 10 and 10°) give comparable results. 

* 

Additional studies showed that the C response corresponding to the 
gain matrix for the feedback configuration / 


T 

y 


(n , q , a j ffl ) 
zcg’ ^cg’ v 5 g' 


was even worse than open loop, having an overshoot of 1. 36. 
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Sensitivity Evaluation 

Because the parameter defining the aircraft equation of motion will be 
changing as a function of mach number and altitude, it is important to deter- 
mine the regions over the flight envelope for which a set of gains will give 
acceptable performance. Also of interest is any modification that can be 
made to the design procedure to desensitize the performance to flight condi- 
tion changes. To illustrate these principles, the feedback configuration 
T 

y =(n , q , a , o> ) was used for designing (from open loop) gain matrices 

zcg eg V g 

for various process representations. These gains were then evaluated using 

data for the cruise, climb, and landing flight conditions. 

For comparison purposes, Table 3. shows the effects of gains designed 

for the cruise conditions evaluated at all three flight conditions. In order 

to attempt improvement of the behavior especially at the landing condition, 

various combinations of process noise (co in eq. l) and measurement noise 

P 

(Y in eq. 4) were introduced into the system model. The corresponding co- 
variances were selected to be proportional to the estimated uncertainty in 
each equation as follows: 

(“ P i 2) “V A p 2(i ’ 3 > £(x p 2 ( j» 


where AA was computed as the average absolute deviation in A over the 
P ~ P 

flight conditions and £L(x 2 (j )) was obtained as the corresponding steady 

XT 

state covariance in x ( 3 ) from a typical run. Note however from Table 10b 

P 

that addition of the process noise computed according to this method resulted 
in a degradation. In fact when combined with measurement noise, the resulting 
gains when applied to the landing configuration resulted in a pair of unstable 
eigenvalues (+.00201 +3 .144). Reductions in the modelled noise variances 
(by 10 to 10 5 ) and retention of only the measurement noise made improve- 
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ments to Table 9 but nothing significantly better than the results in 
Table 9* Similarly various designs using the average process matrices 
for the 3 flight conditions did not appear to be any better than the design 
based upon the cruise condition without any noise. 

Thus as with full state feedback design, the development of a controller 
which takes into account flight condition changes is somewhat of an ad hoc 
process, and further procedures such as those in ref. 17 should probably be 
considered. 
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5 . Conclusions and Recommendations 

Based upon the results presented in section 4 , it can be concluded 

that acceptable flight controllers can indeed be designed using the developed 

reduced state feedback, design program. Typical objectives might include the 

reduction of vertical accelerations and structural loads due to a gust input 

* 

and the response to a C command. 

For a given feedback configuration it is recommended that the design 
process consist of the optimization of a series of feedback configurations 
starting with only one feedback output and progressing one output at a time 
until the desired structure is achieved. Measurement noise should 
be modelled since if present but unaccounted for, severe accelerations can 
result . 


Performance with reduced state feedback controllers can be comparable 
to that achievable by full state feedback systems which in reality cannot be 
designed because of limitations in feeding back the bending effects and in 
designing feedback around the actuators. 

Shortcomings include the excessive computer time requirements (« 1 hour 
for 21 gains and a 17th order system), the existence of multiple minimum 
points, and the sensitivity to flight condition changes. 

Recommendations for future consideration include: 

• Develop parallel (eg. see ref. 18) rather than serial type computa- 
tional algorithms, and implement the design on say the STAR computer 


system. 


2 

Incorporate sensitivity penalty terms of the form (cK . . ) in the 

px oaxj 

performance index, and consider other desensitization procedures as 


per ref. 17. 


• Compare results with those corresponding to a full state feedback 
design implemented with either a linear observer or a Kalman filter. 
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A.l 
A. 2 


where 


where 


where 


Appendix A Program Description 


Program Name : SIRSFB 


Problem solved: 

Process: x=Ax+Bu+G2£+v (A.l) 

ppppppnp 


Reference system: x = A x + B v (A. 2) 

r r r r r s 

Disturbance system: x-Ax+Bv (A. 3) 

n n n n n 

x^ = (NXP x l) plant state 

x^ = {NXR x l) reference state 

x = (NXN x l) external disturbance state 
n 

u = (NUP x l) control vector 
P 

w = (NXP x l) zero mean white plant disturbance with covariance W 
P p 

w = (NWR x l) zero mean white -reference excitation noise with 
2 * 

covariance W 
r 

w = (NXN x l) zero mean white disturbance excitation noise with 
n 

covariance W 
n 


Outputs : 


Z = D x + D x 

p pp p pn n 

Z = D x 

r rr r 

Z^ - (NZP x l) plant output 

Z^=(NZR x l) reference output (NZR = NZP) 


Control: 


u=K y+K y+K y 
P yp p yr •'r yn n 


(A. 4) 
(A- 5) 
(A.6) 


(A. 7) 


y * C x 

P PP P 

y ~ C x 

r rr r 


+ C x + v 
pn n 


n 


= C 

nn 


x 

n 


(A. 8) 
(A. 9) 
(A. 10) 
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y - (NYP x -l) plant feedback vector 
1? 

y = (HYE x l) reference feedforward vector 

y = (NYU x l) disturbance feedforward vector 
n 

J » (HYP x X) zero mean white sensor noise 


Index: J = Lim 
t f -x» 


£[ f 

J o 


[(z - z ) T Q(z - z ) + u T 

p r P ' r p 


B u 3 dt 
P 


'(A 


A . 3 Program Limits 
Variable 


n 


u 


n 


n 


Maximum Dimensions 
18 
12 
6 
6 

12 

6 

12 

12 

6 


A. 4 Theory 

Optimal gains are determined by using the IMSL subroutine ZXPOWL 

(3,4) 

which incorporates the Zangwill-Powell Search procedure. This algorithm 

which does not require gradient computation is such that if the performance 
index were quadratic in the gains, then the search would proceed along a set 
of conjugate directions . 

Following the reading in of the problem description the following 
augmented system is formed: 



■ 

x = A x + B u + v 
P 

y = C x + n 
z = D x 
u = K y 

where 
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(A. 12) 
(A. 13) 
(A. lit) 
(A. 15) 


(A. 16) 


(A. 17) 



(A. 18) 


(A, 19) 


(A. 20) 


(A. 21) 


As shown in Sec 2, the corresponding objective can he written as: 

J = Trace [ (P T Q D + C T K T R K C) P] 

A 


(A. 22) 
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where 



( Q -Q ■ 

\-Q Q y 


K = 


(K K 
yp yr 


K 

yn 


) 


and P satisfies the steady state covariance equation 


(A. 13) 


(A. 24) 


0 = (A - B K C) P + P(A - B K C) T + B K N IC T B T + V (A. 25) 

where R = n) (A. 26) 

V = €(y V T ) 

Thus given initial or perturbed values for K, the- objective J can be 

evaluated by first solving eq. A. 25 for P and then using this value in eq. A. 22.' 

Solution of eq. A. 25 is accomplished using the method of Bartels^ and Stewart. ^9) 

For initialization it is necessary to use a gain matrix K such that all 

eigenvalues -of (A - B K C) are less than zero. If in a subsequent iteration, 

the gain is -perturbed such that one of the eigenvalues is riot negative, then 

50 

the objective is arbitrarily set equal to 10 , thus forcing the optimization 
procedure to backtrack. 

A . 5 Input Format 

Data input consists of system dimensions , defining matrices , and 
various control parameters. If there is no reference system then RXR should 
be read as "zero" and the remaining associated data (AK, BR, CRR, DRR, VR) 
eliminated. Similarly if RXR “0, Then no data cards should be included for 


AN, BR, CPN, DPR, CM, VR. 


All matrices are stored in vector format. 

Rote that the initial gain matrix K must be such that (A - BK) 
is stable. Card format and content are shown in Table A. 1. 
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A . 6 Output 

The following data is printed out: 

© All input data 

a Initial gains and the corresponding eigenvalues, 
a Intermediate values for the objective function. 

• Optimal results (preceded by "OPTIMAL INFINITE TIME SOLUTION") 
consisting of final gains, corresponding eigenvalues, steady- 
state covariance of the state x, and steady state covariance of 
the output z. 

If the initial set of gains is such that the closed loop system is 
unstable, the program will attempt to find a new set of gains as discussed 
in Section 3 * 



Major Subroutines Used 


KEADIN: 

SETUP: 


FUNCTION 


ATXPXA 

SYMSLV 

HSHLDR 

BCKMLT 

SCEUR 

SYSSLV 


EIGRF: 

ZXPOWL: 


Reads in problem data, prints out problem data 
Places the problem into the format described by 
eqs . A . 12 - A . 15 

FZX(XKO); Computes the objective (FZX) given the gain 
matrix (XKO) 

( 19 ) 

Used in the Stevart-Bartels ' ' 

solution of eq. A. 25 


IMSL routine used to compute eigenvalues 

IMSL routine used to perform the optimization 
(described on next page) 


Usage restricted as per letter following ZXPCWL description 
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SUBROUTINE ZXPOWL <F»EPS*N*X*FMIN* ITMAX»tfA» IER). 

I 

-ZXPOWL S/D ’-LIBRARY ,1 

FUNCTION 


USAGE 

PARAMETERS 


F 

EPS 

N 

X 


FMIN 

ITMAX 


WA 

IER 


PRFCISION 

REOtO XMSL ROUTINES 
LANGUAGE 


POWELL'S ALGORITHM TO FIND A (LOCAL) MINIMUM 
OF A REAL FUNCTION OF N REAL VARIABLES 
CALL ZXPQwi. (F,EPS*N»X,FMIN* I TMAX, WA ♦ IEp) 

A FUNCTION SUBPROGRAM WRITTEN 8Y THE USER 
CONVERGENCE CRITERION - SEE ELEMENT 
documentation 

LENGTH OF THE VECTOR ARRAY X (INPUT) 

A VECTOR ARRAY OF LENGTH N. ON INPUT i X IS AN 
INITIAL GUESS FOR THE MINIMUM. ON OUTPUT 
X IS THE COMPUTED MINIMUM POINT 
F (X) - FUNCTION F EVALUATED AT X (OUTPUT) 

ON INPUT = THE MAXIMUM ALLOWABLE NUMBER OF 
ITERATIONS PFR ROOT AND ON OUTPUT = THE 
NUMBER OF ITERATIONS USED 
A VECTOR WORK AREA OF LENGTH N»(N+4> 

ERROR PARAMETER (OUTPUT) 

TERMINAL ERROR = 126 + N ‘ 


ZXPODOIO 

ZXPO0020 

ZXP00030 

ZXPO9O40 

ZXPQ0050 

ZXP00060 

ZXP00070 

ZXP00080 

ZXPO0090 

ZXPOOIOO 

ZXPOOllO 

ZXP00120 

ZXP00130 

ZXP0Q140 

ZXP00150 

ZXPO0X60 

ZXP00170 

ZXPG0I8Q 

ZXP00190 

ZXP00200 

ZXPO0210 

ZXP00220 

ZXP00230 

ZXPO02*»O 

ZXP00250 


N = 1 NO FINITE MINIMUM OBTAINED 
N = 2 F IS LFVEL ALONG A LINE THROUGH X 
N = 4 FaIlURF TO CONVERGE IN ITMaX 
ITERATIONS 

N = d GRADIENT 'LARGE' AT CALCULATED MINIMUMZXP.O0260 
SINGLE/DOUBLE ZXP00270 

UERTST ' ZXP00280 

FORTRAN ZXP00290 

ZXP00300 


CALL ZXPOWL(F,EPS,N,X,imN, ITMAX, WA, IER) 
Purpose 


This routine uses Powell’s algorithm to find a local minimum of a real function of N real vari- 
ables . 

Algorithm 

Let F(X^,X 2 >. . .»Xy) be a function of N real variables X^jX^,. . . ,X^. ZXPOWL seeks a point X*=‘ 

T 

which furnishes a local minimum to the function F at X*, i.e. , 

F (X*) =min (P (X) ) . X in S where S is an open set in E N . 

Note that there are no side constraints so that the problem is simply an unconstrained minimiza- 
tion* 

ZXPOWL uses Zangwill's modification of Powell's conjugate direction algorithm to perform the 
minimization. The algorithm has the notable feature that it will minimize a quadratic form in 
a finite number of steps. 

See references: Zangwill, W, , (1967) “Minimizing a function without calculating derivatives 1 *. 

Computer Journal , Vol. 10, pp. 293-296. 

Powell, (1964) "An efficient method for finding the minimum of a function of several var- 

iables without calculating derivatives". Computer Journal , Vol. 7, pp. 155-162. 

Programming Notes 

The user must furnish the function F as an EXTERNAL FUNCTION subprogram F(X) where X is an N-vec- 
tor of coordinate abscissa. The user must not alter the values in X. 


ZXPOWL-1 



Convergence of the algorithm is defined as 


( 1 ) 




HAXU.,lF(X m ~ 1 )|] 


less than EPS 


If an X 01 is found within ITMAX iterations satisfying (1) X 1 * is accepted as the problem solution, 
and returned to the user. Also the approximate minimum of F, i.e., FCX* 1 ), is returned in the 
output parameter FMXN* Here X 21 denotes the value of X at the m-th iteration. 


Accuracy 

let F(X)=SIN(X)+COS(X) and G(Y)=3Y 3 +5Y 2 +Y+4. Setting X Q =3, Y Q =1, and EPS=10~ 6 , ZXPOWL obtains 
values given in the following table: 


Number of 


Function 

Exact Solution 

ZXPOWL' 

Iterations 

F 

3.926991 

3.927358 

2 

G 

-0.1111111 

-0.1111725 

4 


Note that (1) does not guarantee X? 2 to be close to X*/ It does imply for "well-behaved 11 func- 
tions that F(^) is close to F(X*). 


ZXPOWL— 2 




March 16, 1976 



Professor Howard Kaufman 
Associate Professor 
Electrical and Systems 
Engineering Department 
Rensselaer Polytechnic Institute 
Troy, NY 12181 

Dear Professor Kaufman 

Thank you for your letter of March 10, 1976. IMSL grants 
permission for Rensselaer Polytechnic to use IMSL subroutines 
EIGRE, ZXPOWL, and ZXMIN as part of an application package 
being developed for NASA-Langley. IMSL requests, if possible, 
that Rensselaer Polytechnic only make the object code avail- 
able to NASA— Langley. We also ask that NASA-Langley be _ 
informed that these routines are proprietary and may only 
be used as part of the application package for which they 
were developed. 


Please let us know if we may be of any further service. 


Best regards 





L. L. Williams 
Director, Operations 


maa 


INTERNATIONAL MATHEMATICAL & STATISTICAL LIBRARIES, INC. 


SIXTH FLOOR — GNB SLOG 7500 BELLAIRE HOUSTON TEXAS *77036 [713) 772H927 


Table AM 


45 .” 


Note: Cards for AR, BR, CRR, DRR, VR can be eliminated if NXR=0. 

Cards for AH, BN, CPN, DPN, CNN, VN, GP can be eliminated if NXN = 0. 
Cards for VP should be eliminated if NVP = 0. 


Input Format for Program SIRSFB 


Card 

Column 



Number 

Number 

Description 

Format 

1 

10 

NXP = Dimension of x 

13 

8110 


20 

sr 

NXR = Dimension of x 

- 



r 



30 | 

NXN = Dimension of x 




n 



4o 

NVP = Dimension of v 




p 



50 

NVR = Dimension of v 




r 



60 

NVN = Dimension of v 




n 



70 

NUP = Dimension of u 




P 

I 

2 

10 

NYP =* Dimension of y 

I 8110 



P 



20 

NYR = Dimension of y 




r 



30 

NYN = Dimension of y 




n j 


1 

40 

NZP. = Dimension of z ! 




P 



50 

NZR = Dimension of z 




r 



60 

ITMAX = Maximum No. of iterations 

1 


3 + i 

1 

AP (i, 1) eq. A.l 

8El0.li 


11 

AP (i, 2) 



etc. 

« * * j 


3 + UXP 

1 

AP (NXP, 1) 

8E10.4 


11 

AP (NXP, 2) 



etc . 

« « • 


2 + NXP + i 

1 

BP (i, l) eq. A.l 

8E10.4 


11 

BP (i, 2) 



etc. | 

• * * 


3 + 2 NXP 

i 

1 3 

BP (NXP, l) 

8E10.4 


* I 

BP (NXP, 2) 



1 

etc. j 

i * • • 























Page 2 of Table A.l 


Card 

Humber 

3 + 2 NXP 


Column 

Number 


3+2 NXP 
+ NYP 


3 + 2 NXP 
+ NYP + i 



PP (i, 1) 
PP (i, 2) 

« • • 


CEP (NYP, 1) 
CPP (NYP, 2) 


DPP (i, 1) 
DPP (i, 2) 


3 + 2 NXP 
+ NYP + NZP 


j DPP (NZP, 1) 
DPP (NZP, 2) 


k + 2 NXP 
+ NYP + NZP = 



t NP + NXR 


j AR (NXR, 1) 
AIR (NXR, 2) 


t NP + NXR + i 


j BR (i, 1) 
BR (i, 2) 


• NP + 2 NXR 



BR (NXR, 1) 
BR (NXR, 2) 


Only if TO ^ 0 
Only if NXR y 0 


















r Only if NXR r 0 


Only if NXN f 0 




















Page 1 4 of Table A.I 


48 . 



*T*V 


Only if WXN ± 0 
Only if NYU $ 0 
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Page 5 of Table A . 1 


Card. 

Number 


Column 

Number 


Description 


"Format 


NPRN + i 


1 

11 
etc . 


R (i, 1) 
R (i, 2) 


eq. A. 11 


8E10.U 


NPRN + NUP 


1 

11 

etc . 


R (NUP, l) 
R (NUP, 2) 


8E10.4 


NPRN + NUP 
+ i 


1 

11 

etc. 


Q (i, 1) 

Q (i, 2) 


eq.- A. 11 


8E10.4 


1 
11 
etc . 


j Q (NZP, 1) 
Q (NZP, 2) 


1 
11 
etc . 


j W (i, 1) 
I WW (i, 2) 


CoV of y in eq.A.8 


8E10.4 


8E10 « 4 


NPRN + NUP 
+ NZP + NYP 


1 

11 

etc. 


J W (NYP, 1) 
} mi (NYP, 2) 


1 8E10.4 


5 


NPRN + NUP 
+ NZP 
+ NYP + i 


1 

11 
etc . 


XKO (i, 1} 
XKO (i, 2) 


Initial value of K 
in eq. A. 24 


8E10.4 


















Appendix B 


TIFS DATA (from ref. l4) 

* * d • 

(u, ©, a, 0, n £ , n , n , n^, n^, n ? 

(5 , 5 , 6 , 0 )) 

v sa 5 z* e 5 g' 

^zp’ n zcg i n zt’ n zo)t’ n zFSF, n zast J n zRHT } q cg 


See pages 21, 22 for definitions 



LANDING 


COLUMNS 1 THRU 10 


1 *5 75 18D-01 
-9.483260-01 
-2 .765 220—0? 
0 .0 
0.0 

l .652500+02 

0.0 


-0. 330210+00 
6.770150+01 
-1.300330+00 
1. 000000+00 
0.0 

-4.469700+02 

0.0 


9. 932750-03 
-3.333600-01 
■4. 650270 “0 3 
0.0 
0.0 

6. 509450+01 

0.0 

■1.928430 +0 1 

0.0 

1.051260+01 
0. 0 

4.076500+00 


-6 .721 460 + 01 

0 .0 

6.6 22 9 5 0 + 01 

0.0 

4.41861D + 01 


-4* 125080 + 02 

0.0 

8. 35 66 2D + 0 2 

0.0 

1 . 193480+03 

0.0 

-6.41572D+01 


0.0 0.0 
1.315640+01 -3 .621 300+01 

COLUMNS 11 THRU 14 


-1.40 58 50-02 
6. 220000-0 1 
3. GO 5 5 2D— 0 2 
0.0 
0.0 

-9.374000 +01 
0. 0 

3.77405D+0 1 

0.0 

-2. 361230 + 01 

0.0 

-2* 173790+03 

0.0 

1. 207050+01 


-I .939 840-05 
? .432 410-03 
2.435 070-03 
0 .0 
0 .0 

2 .4 17 73 D >-00 

0.0 

1 .144240+00 

0.0 

-3.9 53 290+00 
1.000000+00 
-I .0 10020+01 
0 .0 

-4. 674 070-01 


1 ♦ 4070 20 —0 2 
-2.712570-01 
2. 411960 -02 
0.0 
0.0 

1.25416D+02 

0.0 

-5.003200+01 

0.0 

1.046400+02 

0.0 

1.64210D+02 

0.0 

-3.761990+03 


* r~ t , s- * 

% - F % + G u 

cj rr AX + 3 LL 


— 1 • 66949D — 0 3 
3.09GUD-01 
1 .48 6 3 70 -0 2 
0.0 
0.0 

-4.7063 20+02 

0.0 

1 . 77 44 40 +0 1 

0.0 

-1.402020+01 

0.0 


-1 .9 19G2D-03 
1 .641650-02 
1 .391 170-04 
0.0 

1 .0 00000+00 
-6 .G4075D+ 00 
0 .0 

2 -08577D+00 

0 .0 

-1 .136 070+00 

0 .0 


-9.732730+00 
-1 .20088D+00 
5.60551 0-03 
0.0 
0.0 

6.795980+ 00 

0.0 

1 .086630+ 00 

0.0 

-2. 016000+00 

0.0 

-4.433250+00 

0.0 

-0.179480-01 


4 . 04 25 ID +00 

0.0 

8.268770+00 


1 .368 91 D+ 00 

0.0 

2 .395 60D+00 


1. 576140-02 
-6. 5966 70-01 
-6. 913040-02 
0.0 
0.0 

9. 175490+01 
0. 0 

-1. 0309 00+03 

0.0 

6. 2198 70+01 

0.0 

4. 666830+01 
0. 0 

-U97382D + 01 


9.2 09560-04 
-7.534060-03 
-2 -103750-03 
0.0 
0.0 

1.525950+00 
1.000000+00 
-3. 5 88 970+ 00 


3. 30895D-02 
1.6780 ‘>0-01 
1. 837090-01 
0.0 
0.0 

1.373110+02 

0.0 

6. 669710+01 


-1 .4 94 570-04 

1 .2 30 92 0-02 

2 .2 2356 0-03 

0.0 

0.0 

-1 .71880 0+00 

0.0 

2.0 2753 0+00 
1 .0 0000 0+00 
-5 .4 2779 0+00 
0 . 0 - 

-2.8 05 89 0+00 
0.0 

1.4 5467 0+00 


0.0 

2.49482 D+00 
0.0 

2*146640+00 

0.0 

-1 .3 76 70 D+00 


0.0 

•1.810250+03 

0.0 

2. 111250+02 

0.0 

6- 2660 9D+00 


3.048490-04 
-9.473790-04 
2 .9 86 22 0—04 
0.0 
0.0 

2.4 6718 0+00 

0.0 

-1.539370+00 

0.0 

1.44632 0+00 

0.0 

1.596200-01 
1.00000 D+00 
-6.526590+ 00 


g' 


7. 28 1000-03 
— 1. 74 570D-0 1 
-9. 26 3920-0 3 
0.0 
0. 0 

6. 1907 6D+01 

0.0 

—3. 06 521D+0 1 

0.0 

1. 804120+01 
0. 0 

-5. 171370 + 00 

0.0 

—3. 58 557D+0 1 


-5.490 120-03 
-1.141 290-01 
-3 .449 84D— 03 
0.0 
0 .0 

4 .485 80D + 01 

0 .0 

-5.776380+00 

0.0 

-2.064 830+01 

0.0 

—9 .960 160+00 

0.0 

1 .927 570+01 


2 .532310— 0 2 
-1.085010-01 
-4.645680-02 
0.0 
0.0 

-1.257650+01 

0.0 

-3.051600+01 

0.0 

5. 188450+01 

0.0 

5.814850+01 

0.0 

1.027340+00 


0.197530-02 

-1.006670+00 

-3.674380-02 

0.0 

0.0 

1 .658000+02 

0.0 

-6.720540+01 

0.0 

6.6 18070+01 

0.0 

4. *02310+01 

0.0 

-3.627090+01 


vn 

H 



LANDING (Cont.) 


COLUMNS 1 THRU 10 


■3# 36 266D— 0 3 
■U 29 77^0-0 2 
0* 68ft* 60 -0 5 
6. 77 89 20 -O 1 
i: S3 7940 -0 1 
1,652070-01 
2.842850-01 
8.881810-04 
3. 48 9110— 0 4 


-3.246 140-02 
-2 .567 800-02 
-3,896 390-02 
-1.824 15D+00 
-4.005 500-01 
-4 .3 12 860-01 
-2 . 0 84 130“* 00 
2 .307690-03 
8 .430 140-01 


3.876640-01 
1 . 25603D— 0 1 
-1.725980 **00 
1.2954 50+00 
6.92 6 3 20- 01 
4. 893570— 01 
—4. 1709 10 + 01 
5.72 6 4 60 + 01 
-8.558460+00 


1 .708 90 D— 03 
1.284 25 D— 03 
1. 163910-02 
-6.762630-02 
-1 .604490-02 
— 1.6 ] 5 60 D— 02 
1.452890-01 
1.246150-03 
8. 48871 0-04 


-3.370910-01 -7 
-1. 203200-01 -3 
-3. 23 12 ID -01 -7 
3.31985D+00 ft 
9.85487D-01 1 

1.03 2400+00 1 

-4 . 1401 2D-0 1 -9 
-4.886040-03 -2 
-1.74 1520-03 6 


.2 70 ft 20- 04 
.9 19300-04 
.2370*0-04 
.302880-02 
.226410-02 
.3 11 880-02 
.407740-03 
.305820-03 
.449730-03 


-1.725700+00 
2. 14306D-0 1 
-1. 3624 20+00 
-4. 0835 50+00 
-5.7194 00-01 
-9, 20 12 50-01 
-5. 776760+00 
5.89 00 30-03 
2. 0982 20—0 3 


-2.605900-03 
6.0 44 39 0-04 
—4 .64*660-03 
-3.7 92 430-02 
—3 .2 64020—03 
-4 .2471] D— 03 
-9.822050-02 
-9.445440-04 
1 .323010-02 


5.720640-01 
-1. 10665D+00 
1. 49 37 80+00 
3.22142D + 00 
2. 42 2410-01 
-7.9020 20-01 
2. 76953D+0 1 
-5.021050-03 
-3.60 5900-0 3 


A = 


COLUMNS 11 THRU 14 


9. 43 2690-0 1 
7. 13556D-0 1 

2. 59 84 50 +0 0 
-4. 66810D+00 
3.065950-01 
9. 62 5770-01 
4. 909330+01 
1. 410880-02 
-7. 359260-04 


-1 .2 58 650-03 
8.933 15D-04 
-2.52? 900-03 
-3 .198 140-02 
-2.213 800-03 
-1 .402 96 0-03 
2 .873 910-01 
-1. 484110-03 
3.131010-03 


1.991650+00 
-2.607760+00 
3.077020+00 
—3. 53070D+0 1 
6 . 404 1 ID +00 
5. 5530 50 +00 
-L. 660240+01 
2.344140-03 
5.94 9 4 90 -04 


1 .948080-03 
—2.6 05 23D— 03 
4*62 2 24 D— 03 
-8.5 1 050D— 02 
5.676970-03 
4.312190-03 
-4.3 52 09 0-02 
1. 740500-02 
-3.^90 rBD-03 


B = 


5.03505D-05 
-3.210950 -03 
-6. 47 4590 —0 3 
-9.299840-01 
-9.41762D-02 
-1.072400-0 1 
*3 00 7D -0 1 
9.541830-04 

3. 343760-04 


1 .12642D T 02 
6.386 120-03 
3 . 099 920- 03 
-1 .2 02 990-01 
-1 .338250-01 
-1 .462070-01 
5 .61460D-O1 

4. B47 980-04 
1 .475 74D—04 


-1.693910-02 
1.068100-02 
-8 . 335 930—0 2 
-1.598350-02 
-1. 154600-02 
-1.00247D-02 
-2.152660+00 
-5.63 24 10-05 
1.690950 -04 


—2.8 0271 D-02 
-3. 109130-02 
-5 .54*970-02 
-1 .8353*0+00 
-4. 153360-01 
-4.389640-01 
-2 .0*634 D+ 00 
2.367020-03 
8.43040 0-01 


7.8 78390-04 
-1.4 7710 0-03 

4.2 80 750-03 

3.9 68030-02 
4. 2 74 500-03 
3.17736 0-03 
1.4 57110-01 
4 .300810-02 

-2 .3 4998 0-03 


vn 

rv> 



CLIMB 

F'% •/- G'u. 
Ax -A So. 


X - 


if =• 

COLUMNS 1 THRU 10 


4 . 8394 60-03 

1 . 1 C 872 D -01 

- 6 . 643860+00 

- 9 . 773390+00 

— 1 . 64 39 80 -03 

- 2 . 96 0260-01 

-1 .001 130 + 00 

8 . 264340 + 0 ! 

- 7 . 991810-01 

2 * 07 57 TO —0 1 

- 3 * 26716 D -0 3 

-2 *094 950-02 

- 1 . 204410+00 

2.90302 D — 03 

1. 1020 3 D -02 

0.0 

0 .0 

1 . 000000+00 

0.0 

0.0 

0 . 0 

0.0 

0.0 

0.0 

0*0 

5. 4999 10 + 01 

1 .766 080+02 

- 5 . 057390+02 

5 . 327280+00 

- 4.7087 60+02 

D . 0 

0 .0 

0.0 

0.0 

0*0 

- 1 . 5 r . 2 ? 4 D -.0 1 

-7 .2 07 210+01 

- 5.505 970 +02 

6.3 22 44 0-01 

1*63 26 20+0 1 

0.0 

0 .0 

0.0 

0.0 

0.0 

6. 76 4 0 40 +00 

7.147 790 + 0 ! 

9 . 070680+0 2 

- 1.2 312211 * 00 

- 1 . 14 06 90*0 1 

0 . 0 

0.0 

0.0 

0.0 

0.0 

2 . 35 01 30*00 

4 *739 560 + 01 

1 . 302060+03 

- 2 . 69096 D * 00 

7.2722 00+00 

Ci.o 

0.0 

0.0 

0.0 

0.0 

- 1 . 087960+01 

-J .940 230 + 01 

- 6 . 537590+01 

- 5 . 17440 D -01 

8.9265 m +00 

COLUMNS 11 THRU 14 




- 1 . 239910-02 

3 . 4 74 600-05 

1 . 341200-02 

2 . 340 P 8 D -04 


6 . 84 5080—0 1 

2 . 533620-03 

- 2 . 909510-01 

- 6 . 902790-04 


5 . 297400-02 

2 .610 340-03 

2 . 292260-02 

3 * 3 44471 J - 04 


0.0 

0 .0 

0.0 

0.0 


0 . 0 

0 .0 

0.0 

0.0 


- 9 . 4244 10 -<01 

2.6 75 080*00 

1 . 338750+02 

2 . * 97600 + 00 


0 . 0 

0 .0 

0*0 

0.0 


4 . 9244 80 +01 

1 .2 27560 + 00 

- 5 . 276470+01 

- 1 . 620800+00 


0 . 0 

0.0 

0.0 

0*0 


- 4 . 604080+01 

~4 .2 80 420 + 00 

9 . 851990+01 

1.5 1 G 04 P + 00 


0 . 0 

1 .0 00 000+00 

0.0 

0.0 


- 2 . 212250+03 

-1 .0 tb 680 + 0 ) 

1 . 421420+02 

1 . 211320-01 


0 . 0 

0 -0 

0.0 

1 . 000000+00 


1. 63 6560 + 0 1 

- 5.129 630-01 

- 3 . 77494 D +03 

- 6 . / 79 46 f ;+ 00 



-1 .5 03 590-03 

1.44 9610-02 

6.9 93350-04 

3 . 1964 90-02 

-1 . 006460-04 

1 .7 28 50 D -02 

- 6 . 914050-01 

- 7 . 813270-03 

9 . 641040-02 

L *3 51 250—02 

1 . 070930-04 

- 6 . 961190-02 

- 2 . 212720-03 

1 . 03 ZZ ID — 01 

2.2 99 33 D — 03 

0 .0 

0 . 0 

0.0 

0.0 

0.0 

1 .000 000+00 

0.0 

0.0 

0.0 

0.0 

- 7.2 36 640 + 00 

9 . 689910+01 

1.563 140+00 

1 . 5 B 60 TO +0 2 

-1 .914 140 + 00 

0 .0 

0 * 0 

1 . 00000 D+ on 

0.0 

0 .0 

2 .215 230+00 

- 1.04 28 OD +03 

- 3 . 7 J 373 D +00 

6 . 344740+01 

2 . 162720+00 

0 .0 

0 . 0 

0.0 

0 . 0 

1.0 00 00 D + 00 

-1 .1 93 200 * 00 

6 . 62 99 3 D +0 1 

2 . 665590+00 

- I . 01 50 80+03 

— 5 .6 J iO 120 + 00 

0 .0 

0.0 

0.0 

0.0 

0.0 

1 .5 20 540+00 

4 . 919790 + 0 1 

2 . 306530+00 

- 2 . 260070+02 

-3 .0 39 86 0+00 

0.0 

0.0 

0.0 

0.0 

0.0 

2 - 502930+00 

- 2 . 308280 + 0 1 

- 1 . 471930 + 00 

— 0. 59 92 S 0 + 00 

1 .5 L 194 0+00 



5*3031 50-03 

1 . 460260-02 

- 3 . 601450-02 

6 . 20502 D -02 

— 2 . 21 56 30—0 1 

-2 . 569350-01 

- 1 . 708990-01 

-1 . 001440+00 

- 1 . 138810-02 

3 . 4 18 750-03 

— 6 . 155750-02 

- 3 . 906300-02 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

7.92 693 D +0 1 

5 . 019 300+01 

- 1 . 627110+01 

1 . 773570+02 

0.0 

0.0 

0.0 

0.0 

- 3 . 922170+01 

- 7.0 78 990+00 

- 3 . 921350+01 

- 7.205 30 D+Ol 

0.0 

0.0 

0.0 

0.0 

2.283800 + 01 

-2 .7 27 68 0+01 

6.6 0 4 4 50 + 01 

7 . 140650+01 

0.0 

0 .0 

0.0 

0.0 

- 7 . 17 6720+00 

- 1 .3 50 64 0+01 

7 . 499560+01 

4 . 717040*01 

0 . 0 

0.0 

0.0 

0 .C 

- 4 . 638350+01 

2 .503 200 + 01 

1 . 343270+00 

-3 . 940100+01 



CLIMB (Cont.) 


COLUMNS 


X THRU 10 


-4. 1794 80-0 3 
-1. 381180-02 
“2 * 306620-04 
—5- 622010-0 L 

— 1 . 297270-0 1 

— 1 • 3989 60 -0 1 
—I ♦ 875000-0 l 

7* 27 7220-04 
2-4009 80-04 


-3 * 405 520-02 
-2.6)4 330-02 
-3 *941030-02 
-1 *956 850 + 00 
—4 *345 360-01 
-4 *5 86 010-01 
-2 *239 820+00 
2 .419 780-03 
6 .961 090-01 


2.535220-01 
1.832130-01 
-1.563260+00 
1. 702860+00 
7. 377420-01 
6*479810-01 
-4*510460+01 
5.725610+01 
~7.G69270<00 


9.319780-04 
1 .526750-04 
5 .946100-03 
-4.205230-02 
-1 *027640-0? 
-1. 039530-02 
8.662560-02 
9.325520-04 
5.45871 0-04 


-3.384790-0 1 
— 1.2101 2D— 0 1 
-3.305160-01 
3.320000+00 
9.013700-01 
1.028350+00 
— 5.30558D— 01 
-4.777590-03 
-1.4194 30-03 


-7.153300-04 
-4 .060670-04 
-7 .6 69 960-04 
0 .844 91 0-02 
1 .2 0799D-02 
1.377630-02 
-1. 16773 D-0? 
-2.308880-03 
5 .3 23 770-03 


-1. 726710+00 -2.714530-03 
2. 1 3909D-0 1 6.164 390-04 

-1. 361060+00 -4*701060-03 
-4.9Q134D+00 -3*954950-02 
-5.864380-01 -3*277390-03 
— 9. 434440 -01 -4.20849D-O3 
-5.901390+00 -1. 047020-01 
5. 8671 90—0 3 -9.650270-04 
1* 764450-03 1 .092430-02 


5* 747460-01 
-1. 1060 1D+00 
1*40 64 20+00 
3*008670+00 
1*95 2930-01 
-8.388390-01 
2*801220+01 
-6. 223960-03 
— 3* 166090—03 


COLUMNS 11 THRU 14 


-9.400950-01 
7* 138500-01 
-2* 596340 +00 
-4.619530+00 
3. 195790-01 
9*769890-01 
5. 113160+01 
1* 510010-0? 
-7. 277420-04 


-1.2 57 090-03 
9.173 540-04 
-2 .544290-03 
-3.468410-02 
-2 .517 790-03 
-1*704 330-03 
3.071900-01 
-1 *471400-03 
2. 595 760-03 


1*995270+00 
-2.612250+00 
3.06 79 50+0 0 
-3.552920+01 
6*4054 8D+00 
5*556200+00 
-1*60 54 5D +01 
1.643970-03 
3*458920-04 


1 *983550-03 
-2.71752 0-03 
4.791590-03 
-0.884970-02 
5.83634D-03 
4.43673D-03 
—4 .436080-02 
1 .739400-02 
-2*891880-03 


B = 


-6. 15895D-05 
-4.260210-03 
-7* 58 1120-0 3 
-1* 195340+00 
-1- 195010-01 
-1.362340 -01 
-3*341390-01 
1*203310-03 
3. 53390D— 0 4 


— 4. 443 03 D— 03 -2.20 1610-02 -3.069090-02 
-3.320Q1D-03 1.047950-02 -3.378210-02 

2 .6 17 10D— 03 -1*123260-01 -5.940400-02 
-1*753220-01 -2.259200-02 -1*971300+00 
-1.845530-01 -1.789170-02 -4*437910-01 
— 1 .990 040—01 —1*6063 2D —02 -4*608520-01 
7 .4 70 230—01 -2*779320+00 -2.253930+00 
6 .4 19 7? 0-04 - 5 . 2 7 56 2D -0 5 2 . 4 79 26 0-03 

1 .656370-04 2.054570-04 6.961350-01 


8.320)60-04 
-I *4 62 15 0-03 
4*3 6216 0—03 
4*2 8932 0-02 
4.6 63 300-03 
3*5 7760 0-03 
1 .5 26 820-01 
4*2 97620-02 
-1 *9 48910-03 


VJl 

-P- 



CRUISE 


« = Fx + Gu 


COLUMNS 1 THRU 10 
-1.64.3150-02 -L .0 25 620-03 

1 , 4669 QD +00 

y 

-9.805670+00 

~ f\X -h 

-1.35 7130 -04 

3 tM 

—4 .0 78 260-04 

1.612190-03 

1.2 46 260-04 

2. 03 8930-02 

6.9 14490-05 

-I. 51308D-0 1 

-1 .890 560 + 00 

1.487160+02 

9.724760-02 

3.705740-0 1 

3 .132 870-02 

-1.4295 30+00 

-1 .1 52 760-02 

-3. 3012 5D-01 

2.00295 0-02 

3. 329960-04 

-5.426 990-02 

-2.141110+00 

—2 .070420-04 

“4. 236620 “03 

2 .362 820-04 

— 1 . 3891 7D— 0 1 

-4.0 45 820—03 

3. 50 33 60 -01 

3 .7 89 83 D— 03 

0.0 

0.0 

1 . OOOOOD +0 0 

0.0 

0.0 

0 .0 

0. 0 

0.0 

0.0 

0.0 

0. o 

0.0 

0.0 

n.o 

0.0 

1 .000 000+00 

0.0 

0.0 

0.0 

0.0 

2. 16902D+0 1 

3.291 eiD+O? 

-1. 114090+03 

-7.743670-01 

-5. 155830 + 02 

“1 .229020+01 

2. 361470+02 

2.1 39 40 D+ 00 

4. 249500+02 

-4.8 40810+00 

0.0 

0.0 

0.0 

0.0 

0.0 

0 .0 

0.0 

1 .0 00000+00 

0.0 

0.0 

—5* 77646D+00 

-l .385 740+02 

-9. 960420+02 

-4.6 59160-02 

2. 17Q24D +01 

4 .0 30070 + 00 

-I. 120780+03 

-5.285410+00 

9.9680 10+01 

4.0 55 93 0+00 

0. 0 

o.o 

0.0 

0.0 

0.0 

0 .0 

0. 0 

0.0 

0.0 

1.0 00 00 D+ 00 

-6. 40 4970+0 0 

l .416590+02 

1.743990+03 

1 .4 36560-01 

-4.839390+00 

-2 .264420+00 

1. 5023 60+02 

4. 82396 D+ 09 

-1.975010+0 3 

-U .447 03 0+00 

0. 0 

0.0 

0.0 

0.0 

0.0 

0 .0 

0.0 

0.0 

0.0 

0.0 

-4. 08 6970*0 0 

8 .974080+01 

2. 58232D+0 3 

3.80319D-01 

3. 470860 + 01 

2 .6 84 77D+ 00 

1. 140330+02 

4 .4 57200+00 

-5. 33 34 50 +0 2 

—4 .9 1097D+ QO 

0.0 

0.0 

0-0 

0.0 

0.0 

0 .0 

0.0 

0.0 

0.0 

0.0 

—4. 409770+0 0 

-6 .139310+01 

-7. 350040+01 

1 .022680-01 

2. 1792 90 +01 

4 .687640+00 

-6.939200+01 

-2.446970+00 

-3.4996 70+01 

3-0 93670+00 


COLUMNS 11 THRU 14 


4.040590-03 
1 ♦ 644920+00 
2*13 12 5D— 0 1 
0* 0 
0*0 

-1.80990D+02 

0*0 

1.7 f 306D +0 2 
0*0 

-2.73 3420 *02 

0.0 

-2.495750+03 

0.0 

5. 901700+01 


1 .908 190-04 
2. 040 020-03 
3. 126 16D-03 
0.0 
0 .0 

5.5 98 53D+00 

0.0 

2.209240+00 

0.0 

-7.700 020+00 
1 *000 COD + 00 
“1 .765 060 + 01 
0.0 

“1.046 090+00 


5. 790610-03 
-7. 03 1 7 1D—0 1 
1.8639TD-02 
0.0 
0.0 

3.079370+02 

0.0 

“1.21 1300 +02 

0.0 

1.43&96D +02 

0.0 

1.084420+02 

0.0 

-3. 862060+03 


2.800990-05 
-1. a 62710-03 
3.144230-04 
0.0 
0.0 

4.210350+00 

0.0 

-2.0 11230+00 

0.0 

2.686830+00 

0.0 

3 .701070-01 
1 . 000000+00 
-9.214080+00 


G 


/ 


1 .646420-03 
-6 .677 140-01 
5.192 630-03 
0.0 
0.0 

1.72340D+0? 

0.0 

-1.560 640+01 

0 .0 


— 1.5110 ID— 02 
-7. 930650—01 
-1. 92 5 7 00 -01 
0.0 
0.0 

-4.386060+01 

0.0 

-1.154600+02 


'1. 154920-02 
-6* 63 97 5D—0 1 
*3. 60 77 ID— 02 
0.0 
0.0 

2. 359920+02 

0.0 

■1. 17499D +02 

0.0 

7.050930+01 

0.0 

•1.72 323D+0 1 

0.0 

1 .376020+0 Z 


-9 . 1 23 45D+0 1 

0.0 

-5.239 270+01 

0.0 

7-373 73D + 01 


0.0 

1.962050+02 

0.0 

2. 167960+02 

0.0 

3. 2376 70 +00 


7 .969130—03 
-1.914030+00 
-7.911760-02 
0.0 
0.0 

3 . 29756D+ 02 

0.0 

-1.3 07240+02 

0.0 

1.41624D+02 

0.0 

6.98583 0+01 

0.0 

-0.151310+01 


\_n 

\J1 



CRUISE (Cont.) 


COLUMNS 1 THRU 10 


2- 98 3130-03 
- 1 . 054150 - 0 ? 
7.91 0930-03 

2. 014950-01 
■7. 231660—02 
■ft. 121460-02 
1.635900-01 
4.913290-04 
8.055670-05 
■6, 613020 +02 
■6. 701370+03 
7.459640+02 


-6-700 040-02 
—5*045 850—02 
—7 *694 300-02 
-3.753970+00 
-0.016420-01 
-0.45X010-01 
—4 * 336 690+00 
6.120840-03 
3 , ft 37 91 D— 01 
3.001090+02 
-1.264 19D+03 
4.086 370+ 03 


1.920440-01 
1 . 4890 8D— 0 1 
— 2.68 54 30 + 00 
4 . 99 9 7 50 +00 
1.473B80 +00 
1. 345750+00 
—8. 79 1 36D +0 1 
5.720030+01 
-3.900330+00 
1.320900+03 
—4. 821820 +02 
-1. 18 50 ID +0 4 


-2 .78703 0-04 
-1 .120240-04 
-7.222440-04 
6.526090-03 
1 .391750-03 
1.427460-03 
-1* 156550-02 

2.858900-04 

9 .9 30 01 0—05 
3 .961 5BD— 01 
2.43285D+00 
—9.5 83470+00 


—3 . 3765 £0 — O 1 
-1.221990-01 
-3. 526530-01 
3.81 6250 +00 
1.057490 +00 
1. 108350+00 
-1.243270+00 
-7. 1846 80-03 
— 1.1686 00 -03 
1 . 4291 80 + 0 4 
1.076800+05 
-2.735650+03 


-6 .670 640-04 
-5 .3 58 27D-04 
-6 .4 55 060-04 

1 .5 50150-01 
2 .169070-02 

2 .320 090-02 
-1 .705 210-02 
-2.429250-03 

Z .9 03 080-03 
4. 342550+01 
2.4 24 270+02 
-4.262280+01 


-1.708850+00 
2. 1083 3D— 0 1 
-1.4364 ID +00 
-6. 744070+00 
-9. 11803D-01 
-1. 2871 40+00 
-ft. 0765 1D+00 
1.010510-02 
1. 9491 10-03 
2. 794220+04 
„ I. 3955 6D +04 
“ 5 . 06 36 00+04 


-2 .501200-03 

9.1 72 69 D— 04 
-6 .8 9204 D— 03 
-5 *0 49 B4C—02 
—4 #2 05 74D— 03 
—5 .5 01 40 D— 03 
-1.9 00 89D-01 
— I . 1 2338 D— 03 
5 .97230D-03 
4.095560+01 
5.254860+01 
5 .4 17 850+01 


5.6993 10-01 
- 1.1070 1 D +00 
1.676630+00 
7 . 15 OZ 10-0 1 
-2.916800-01 
— 1. 33 65 20+00 
3.740090+01 
-7. 09 3130-0 3 
— 3. 314560-03 

1. 5430 60+04 

2. fil 48 50+04 
9.955980+04 


COLUMNS 11 THRU Z4 


—9 « 92 3910-01 

7. 167980-01 
—2* 48 5170+00 
-3.238900+00 
6. 28 8470-01 
1. 304110+00 
6.006890+01 
2.370030-02 
-1,40 5920-03 
5.160620+03 
2. 04366D<04 
-5.731640+04 


-1 .7 07 060-03 
9.086890-04 
-1 .671720-03 
-6.427 81D-0Z 
-4 .1 72 030—03 
-5 .3 54 61 D— 03 
5 * 1 3ft 220-01 
— 1 .301 790-03 
1.482 200-03 
5 .926 570+00 
4.234 890+01 
-3.747 420 + 01 


2.01600D+0C 
-2.610600+00 
3.090960 +00 
-3. 79 60 20 +01 
6. 112000+00 
5. 23658D+00 
-1.607390+01 
1.654990-03 
4.179590-04 
6.322120+04 
2. 316900+05 
5.01225D +03 


2.273750-03 
-3.039000-03 
5.993180-03 
-1.285070-01 
6.107170-03 

4. 2 31 68 D— 03 
-7.373360-02 
1 .734 370-02 
-1.600 560-03 
8.647430+01 
2 .827580+ 02 
-1 .996740+01 



B 


3. 09 OU 20-0 3 
-1. 1973)0 -0 2 
-2.214770-02 
-3. 551740+00 
-3.584000-01 
—4. 0843 2D— 0 1 
-1. 118340+00 

4. 77344D-03 

8. 271010-04 
4. 570190+02 
1. 02 3840+0 3 
2. 12 2840 +0 2 


8 .748 05D-03 
-4 . 504 350-03 

9. JT 47 1 20-03 
-5. 014470-01 
-5.316420-01 
-5 .784420—0 1 
2 * 636 320+00 
2 *4 21 340—03 
3.048790-04 
1 .604960+02 
—2 . 4 78 07D + 03 
2.241750+02 


-7. 733450-02 
2. 877550-03 
-3.676010-01 
-1.334880-01 
-8. 099510-02 
-8. 476100-02 
-8. 142940+00 
1.260000-04 
6. 36364D -04 
2*40 50 90 +03 
9* 47424D +0 3 
-4.023990+02 


-4.292740-02 
-5.172040-02 
-1.090140-01 
-3.764160+0 0 
—8 .049000—01 
-8.506 20D-01 
-4.374540+00 
6.205920-03 
3.833120-01 
—1.173 95 D+ 04 
-1 .227270+05 
1.543960+03 


8.198430-04 
-1 .56534 0-03 

5.9 6560 0-03 
ft .6 02 75 0-02 
1.07766 0-02 

9.9 40160-03 
2 .3 96 85 D— 01 

4.2 6775 0-02 

-1 .12310 0-03 
3.0 1073 0+01 
2 *2 1 592 0+ 01 

4.9 80980+01 



Appendix C 


Publications and Presentations 

Presentations -which incorporate a publication in a preprint volume. 

"Computation of Output Feedback Gains for Linear Stochastic Systems 
Using the Zangwill-Powell Method", I 977 JACC, San Francisco, CA, 
pp. 1576-1581. 

"Application of Stochastic Optimal Reduced State Feedback Gain Com- 
putation Procedures to the Design of Aircraft Gust Alleviation 
Controller s'", IFAC VII World Congress, Helsinki, June 1978. 

"Application of Stochastic Optimal Reduced State Feedback to a 17 th 
Order Aircraft Model," "Optimization Days, 1978," Montreal, May I 978 
(Abstract published). ‘ 


Relative Master Project Reports. 

"Comparison of Stabilizing Subroutines", RPI, ESE Dept., P. Baratta, 

Aug. 1976. 

"Reduced State Feedback Gain Computation Using the Sequential Uncon- 
strained Minimization Technique", RPI, ESE Dept., J. Yip, June 1977* 

"Investigation of the Design and Performance of Reduced State Feed- 
back Controllers", RPI, ESE Dept., K. Sobel, May 1977- 




Table i 


Effect of various controller configurations 
Initial gains set equal to zero, sensor noise not accounted for in design 


Vertical Acceleration 


mas values 


(ft/sec“ ), 



No sensor 
' noise ' 

mt. sensor 
noise 

Feedback 

Configuration 

gain matrix 

eigenvalues 

J* 

n zl n z 2 

n _ n _ 

zl z 2 

Open loop 

[ 0 ] 

t 

-1.45ij .660 
-.00705 ±j .823 
- 20 . 

-4o. 

.0722 

9.75 9.28 

- 

3 

ri 

CD 

n 

G* 

0 . 10 T 0.0997 0.228 

-193 -53-6 -67.9 

- 20.9 ±0 43.3 

- 18.9 

- 2.07 

- 0.0607 ±3 0.0537 
-0.2784 

.0190 

3*40 3.43 

1 

63.8 54.5 

<1,9, a-a 
& 

- 1.82 -O .188 0.0262 

- 13.1 -ia . 3 - 4.02 

- 36.6 

- 11.7 ±3 10.8 
- 0.223 

- 2 . 6 l 

- 0.050 

-0.2784 

.0150 

I .98 1.92 

7.78 5-84 

<1,0, a, a 

- 0.310 - 0.191 0.0442 - 0.219 
- 1.02 -52.5 -4l6 3.14 

- 20.2 ±j 63.7 
-18.4 
-3-91 
- 0.222 
-0.0377 
-0.2784 

.000528 

2.78 2.52 

71.4 60.7 

1 , 0 , a-a, 

5 

a 

s 

-0.575 0.286 0.0175 -0.0427 

-150 - 83.3 -4.25 0.530 

-20.4 ±j 30.4 

-20.9 

- 1.11 

-0.139 

-0.0995 

-0.2784 

.0143 

2.83 2.78 

46.7 39-5 


VJl 

MD 














Table 2 

Effect of various controller configurations 


Vertical acceleration. 


Initial gains set equal to zero, sensor noise accounted for in design 


rms values 

Ho sensor 
noise 


ft/ sec 2 . 

sensor 

noise 


Feedback 

Configuration 

gain matrix 

eigenvalues 

J* 

n n 

zl 

n „ 
z2 


n zl 

n 

z2 

Open loop 

% 

[0] 

- 1.45 ±5 .660 

-.0075 ±J .0823 
-20. 

-4o. 

.0722 

9-75 

9.28 


- 

- 

<1,6, a 

0.00721 - 0.0299 - 0 . 031*3 
- 0.0239 - 0.0204 - 0.539 

- 39.85 

-20.05 

- 1.48 +j .888 
-00360 +j .0877 
-O.278 

.0667 

9.22 

8.86 


9.23 

8.86 

q,0,a-a 

S 

o.oo 4 n -0.0368 0.0119 

-0.0102 0.0175 -1.18 

- 39.7 

-20.0 

- 1.57 +A . 84 6 
-.0444 + j .0822 
- 0.2784 

. 0478 

4.43 

4.25 


4.42 

4.25 

q,6,a,a 

g 

0.00518 -0.0661 -0.186 -0.l4i 
- 0.0435 -0.0201 - 0.140 1.18 

- 40. 0 
-20.1 

-1. 4 o +j 1.22 
-.0575 +j .0913 , 

- 0.2784 

.0503 

4.74 

4.23 


4.75 

4.24 

q,0,a-a , a 
g g 

.00468 -0.0398 -0.00984 -0.0451 

-0.0100 - 0.00753 - 0.930 0.794 

- 39.8 

-20.0 

-1.52 +j .886 

-.0470 +j .0832 
-.2784 

.o4i4 

6.52 

6.39 


6.51 

6.39 


cr\ 

o 














Table 3a Effects of Initial Gain Selection for J 


y = (u, e, n 5 4__, a , w_) 


Initial 

gains 


OPTIMAL 

GAINS 

FOR 

n } q ,a , 
zcg eg v 


OPTIMAL 

GAINS FOR 

n , q. ,a , 
zcg eg’ v’ 

W , u, n . 


Gains 

eigenvalues 

r-1.69 

-6o4 

3.42 

.0355 

2.84 

-.674 

-2. 19±,}65. 4 

-24.2*31.46 







-10.8*350.3 

-.00595*31.23 

.565 

685 

.568 

1.48 

3.71 

1.15 

-1.78±j37.9 

-4l.8 



1.42 




-2.93*332.3 

-2.35 

-.0150 

30 

-.197 

-.776 

.236 

-8.73*320.3 

-.348 

♦105 

-25-9 

.870 

-.235 

1.71 

-.345 

-3.15*363.4 

-.00320*3.00989 







-9.60*349.6 

-40.0 

.106 

38.5 

-2.93 

2.54 

3.93 

.932 

-1.07*336.1 

-12.9 


-1.62 

2.06 




-5.67*332.2 

-.316 

-.791 

-1.30 

.0656 

-.00322 

-5.30*323.8 

-22.2*312.6 




.402 

1.78 

1.17 

.235 

. 11.2 

-3.91 

-2.64*363.9 

-4.22*32.91 


18.9 


2.47 



-10.1*349.3 

-40.5 

.0517 

-3.72 

2.83 

1.45 

-3.36*336.1 

-.00152 

-.0208 

54.8 

.699 


-4.77 


-.791*334.9 

-.1984 

-1.26 

1.79 

-8.31*322.6 








-23.8*38.21 















Table 3b Effects 

of Initial Gain 

Selection for J , 




Covariances 


Initial 

n z 

n 

BMR 

SR 

p 

zcg 



Gains 

(g 2 ) 

(g 2 ) 

(n 2 -m 2 ) 

d d 
n -m 


OPEN 

LOOP 


OPTIMAL 

GAIN'S 

FOR 

n . q 
zcg’ H cg 3 

a , \>r 
v 5 g 


1265E-03 


.1236E-03 


.1126E-03 


. 2019 E -03 

. 7485 E- 1 -IO 

.7171EH-08 

. 4855^00 

•5344E+00 

. 4894E-02 

.2022E-03 

. 7193 E +10 

. 9 M 1 OE +08 

. 4277E-01 

. 1603 E +01 

.3454E-02 

.1760E-03 

. 706 OE+IO 

.9750EH08 

. 2398 E -01 

.1752E-I-01 

. 3454 E -02 


























Table 4a Feedback Design Study for J Q , Cruise Condition 
Q = DIAG (10 6 ,10^,0) R =» DIAG (4. ,6. ,25.) x 10 -6 


OUTPUT 

FEEDBACK 

CONFIGURATION 

J 

if 

Iter- 

ations 

K such that <$SA 

\ z 

6 e 

= -Ky 

EIGENVALUES 







-4.90*361.6 

-I.52i3i.66 







-ll.4±j50.9 

-0.00689±3 .0465 

OPEN LOOP 

7477 


0 



-l.76±j4l.5 

-40. 







-3.021332.7 

-30. 







-5.891322.3 

-20. 




I" " ' 1 

1.047 



-1.301364.4 

-2.301323.4 







-9.501350.4 

-1.2413 3. 9 5 

n 

1132 

9 

-7.698 



-9.081335.6 

-0.0063l3O.0235 

zcg 






-3.461J34.1 

-40.4 




4.889 



-26.Oi35.98 





3.063 

0.3072 


-1.711365.2 

-1.4H33. 61 







-10.7±350.5 

-0.0064l30.0292 

n ,W : 

24i 

13 

-1.670 

2.253 ' 


-2.861338.4 

-4l.0 

zcg g 






-3.921331.5 

-26.5 




2.253 

-0.00839 

1 

-7.OI+320.9 

-22.9 




O.527 

0.230 

i 

-0.306 

-30.1+363.2 

-3.90*324.2 







-9.35*349.6 

-0.00756130.0187 

n ,W ,q. 

169 

22 

-4.02 

2.37 

2.85 

-0.978*336.0 

-39.8 

zcg’ g eg 






-7.001331.8 

-11.6 




2.17 

0 . 0229 

-1.54 

-22.5±313.9 

\ 

-2.53 




OUTPUT . . . 


#Iter , . 


K such that . * * . 


Eigenvalues 







n ,W ,q 5 

ZCg g Cg 

a 

V 

153 

17 

.578 -4.12 0.0648 11.46 
-3.04 3-04 5-41 -1.31 
.00571 0.912 -0.967 -2.25 

-3.46±j63.6 -7.44±j22.1 

-10 . 5±j 49 . 5 -2.82ij3.32 

-2.66ij36.2 -.00853ij0.056l 

-.00895±J34.2 -39-9 

-26.9±j7.77 

n ,W ,q , 

zcg’ g ^eg’ 

-2 

a ,n _xlO 

68.9 

26 

.771 -4.64 .-.539 9.25 -51.9 

-12.5 2.33 .435 3.72 94.2 

-.405 1.47 -.894 -3.37 19.4 

-.00845*,) 56. 6 -.00760±j0i0748 

~10.7±j49.4 -103 

-6. 07± j 34. 6 -23.6 

-2.24ij32.6 -16.1 

~8.o4±j20.89 -2.51 

n zcg J ^g’ q eg 5 

a v’ n zvT xl0 ~ 2 ’ U 

67.7 

9 

.758 -4.67 -.634 9.24 -53.4 0.210 

-12.5 2.31 .406 3.72 94.1 0.010 

1 

1 

-.236 1.49 -.629 -3.37 19.9 0.100 

~.00l62±J56.6 -3.50±J2.64 

-I0.9±j49.6 -106 

-6.17iJ34.9 -22.5 

-2.31iJ32.4 -23.5 

-7. 87+ j 21.1 -0.0202 

-0.00217 

^cg’^g’^cg’ 

a v> n a »T xl52 > e 

65 .k 

31 

1.27 -4.74 -1.56 9.10 -53.5 1600 

-13.5 2.4l 1.93 3.75 100.0 -826 

-.526 1.48 -.925 -3.39 17.6 -231 

-.Il8ij56.8 -106 

-10.7±j49.2 -26.2 

-5.08±J33.8 -0.0000624 

-10.81J21.8 -0.0916 

-6.91±J7.19 -6.68 

-3.29±J32.9 

n zcg’ W g ,<1 cg’ 

V^vT* 1 ^’ 

BMRxlO" 5 

68.1 

ll 

.790 -4.65 -.650 9.24 -53.0 .0111 

-12.6 2.34 .507 3.69 94.2 .00573 

-.294 1.47 -.562 -3.35 20.5 -.00758 

-.000711ij 56.7 -2. 98+ j 3. 07 
-10,9±J49.6 0079 4iJ 0.0734 

-6,09ij35.0 -106 

-2.27iJ32.4 
-7.86±J21.2 

~23.6ijl.62 

























Table 4b Feedback Design Study for J^, Cruise 
Q « DIAG (10 6 , 10 6 , 0) R = DIAG (4., 6., 2 

Condition 
>5.) x 10" 6 


F 

FEEDBACK 

CONFIGURATION 

” 4p (« s ) 

COVARIANCE 

■W^ 

COVARIANCE 

BMR(n-ttl) 2 

COVARIANCE 

S R Cn-m) 2 

COVARIANCE 


<S z (dg) 2 . 

COVARIANCE 

Mag) 2 

© r 

COVARIANCE 

OPEN LOOP 

0.743OE-02 

0.7525E-02 

0.1128E+11 

0.1121E+09 

0 

0 

0 

n 

zcg 

0.8371E-03 

0.1427E-02 

0.4442E+10 

0.4128E+08 

0 . 1808 E -02 

0.9H2E-01 

0.3846E-01 

n , ¥ 

zcg’ g 

0.2072E-03 

0.2750E-03 

0.7433E+10 

0.9555E+08 

0 . 558^-01 

0.1504E+01 

0.2522E-02 

n , ¥ , q 

zcg g ^cg 

0.1239E-03 

0.2146E-03 

0.7003E+10 

0.9455E+08 

0.1915E-01 

0.1579S+01 

0 . 3509 E -02 

n ,¥ ,q ,a 

zcg’ g eg’ V 

0.1402E-03 

0.1649E-03 

0 . 6696^+10 

O. 987 OE+O 8 

0.4223E-02 

0.1911E+01 

0.3829E-02 

n zcg’ W g’^cg’ a v’ 

■w* 10 " 2 

0.6874E-04 

0.6902E-04 

0.4774E+10 

0 .1088E+09 

0.4989E+00 

0.3356E+01 

0.1691E-01 

n zcg’^g’ q -cg’ a v’ 

n zvT xl0 " 2 ’ U 

0.6616E-04 

0 . 6917 E -04 

0 . 5220E+10 

0.1106E+09 

k 0.8124E+00 

0.3422E+01 

0.225.7E-01 

zcg’ g’^cg’ V* 

n zwT xl0 ~ 2 ’ 0 

0.6310E-04 

0.6773E-04 

0.5195E+10 

l 

0 .1084E+09 

0.2625E+01 

0.4477E+01 

0.3379E-01 

n zcg’ W g’ 4 cg’ a v’ 

n zrt x10 ” 2 * 

BMRxlO" 5 

0,666tE“04 

i 

| 

0.6959E-04 

0.4882E+10 

0 . 1098 E +09 

1 

i 

0.4650E+00 

0.3413E+01 

0.1652E-01 


G\ 















































































Table 5a. Feedback Design Study for Cruise Condition 


OUTPUT 

FEEDBACK 

CONFIGURATION 


n ,W 
zcg’ g 


n S W ,q 
zcg g eg 


n »¥ ,q , 

zcg g eg’ 




.564E+l6 


.202E+16 


.214E+15 


.109E+15 


Q = DIAG (0.,0.,10°) 


IC such that 6 SA 


R = DIAG (U.,6.,25.) x 10“ 



“ -Ky 



-5-99 

- 6 . 0 k 

5.92 

-1.44 -6.37 

-1.44 5.56 

- 1.02 .755 


1 .56 -1.37 -3.29 
-7.92 0.307 0.538 
-3-90 0.549 -.434 


-1.64 -5.69 
-5.88 2.16 
-3.24 -.515 


-1.69 23.3 
-2.41 -33.0 
-.960 -5.31 


EIGENVALUES 

-4.9±j6l.6 

-1.53151.66 

-11.4±j50.9 

-.00689+50. 0465 

-17. 615 41.5 

-40 

-30.21532 .7 

-20 

-5.89±j22.3 

-30 

-11.3±557.2 

-24.2±j5.19 

-1.10±j49.2 

-1,20153.83 

-l6.i±j45.9 

-. 00629± j. 0230 

-l.82±J33.4 

-38.4 

-.595*524.6 


-5.8615 59.6 

-.00245±j .0883 

-11.7*551.3 

-39.3 

-1 . 33*5 43 . 2 

-31.2 

-2.66+533.2 

-21.1 

-4.99*523.1 

-2.48 


-0.48 


■l.l6+}64.5 

-.00664+J0.759 

■11.8±j49.4 

-41.3 

-.0347*542.8 

-29.8 

•5.21+532.0 

-19.6 

5-45±j25 .6 

-7.81 


-4. 43ij6l. 4 

-.0082515 .0754 

-11.4+550.8 

-15.6 

-0.720+542.9 

-9.94 

-.728±531.8 

-5.58 

-38.7 


-25.9 


-8.811525.6 






















Table 5a (Continued) 


OUTPUT 

FEEDBACK 

CONFIGURATION 


n zcg jW g’%g* 

VW 10 " 


n ,4 , 

zcg’ g’U:g’ 


n ,¥ ,a , 
zcg g eg’ 


a , d 


n ,w ,q , 
zeg’ g’^cg’ 


BMR x 10 



0 . 600E+14 


K such 

that <5 sa 

6 Z - -Ky 



-1.07 

-5.56 -1.72 

25.1 

-7.82 

-2.69 

2.16 -1.36 

-33.0 

35.8 

-2.87 

.488 -.967 

-6.08 

5.63 

-1.66 

-5.65 -1.49 

23.4 

.00986 

-4.58 

2.16 -2.07 

-32.2 

-.00108 

-3.46 

,486 -1,00 

-5.80 

-.0000129 


EIGENVALUE 


9.31±j59.3 

-.00843*3.0746 

•10. 7*3 49-8 

-48.7 

■.0142+J43.8 

-15.5 

■.8l3±J32.Q 

-8.79 

8, 11+ j 24. 6 


12.OiJ5.l6 





-4.93±j6l.2 

-11.48*350.9 

-.0455^42.8 

-.6l7±J31.8 

-9. 3U±J25*3 

-ll.9-Jl.35 


-5.0l±3oi.lv 

-11.5±J50.8 

-.000188+J42.6 

-.576+331.8 

-9.44+325.4 ■ 

-lO.9t34.93 


-5.31±36o.t 

-10.8+351.2 

-,ooo466±j42.7 

-36.6t32.94 

-.6711332.0 

-10.21J28.3 


-.008o4±j .0787 

- 38.8 

-26.4 

-6.05 









Table J? 1 ? Feedback Design Study for 3. , Cruise Condition 

/• J- ^ 

a - DIAG (0., 0., 10 ) R - DIAG (U.« 6.» 25.) at 1(T' 


FEEDBACK 

CONFIGURATION 

V 8>2 

COVARIANCE 

n (g) 2 

zcg 

COVARIANCE 

BMR(nm) 2 

COVARIANCE 

S R (nm) 2 

COVARIANCE 

4 SA (d8)2 
COVARIANCE 

e(&&) 2 

4 j 

COVARIANCE 

5 (dg) 2 . 

COVARIANCE 

OPEN LOOP 

0.7U30E-02 

0.7525E-02 

0.1128E+11 

0.1121E+09 

0 

0 

0 

n 

zcg 

0 . 3732E-02 

0.1582E-02 

0.1+039E+10 

0 . 1+331E+08 

0 . 5812E-01 

0.5170E-01 

a. 5389E-01 

n , ¥ 

zcg g 

0.18U8E-01 

0.8678E-02 

0 . 1+273E+09 

0.1210E+09 

0.1291E+02 

0.951+8E+01 

0.1734E+00 

n j¥ ,q 

zcg’ g eg 

0.3757E-01 

0.2232E-01 

0.2173E+09 

" ' ' 1 

0.2851+E+08 

0.3537E+01 

0 4 1010E+01 

0.21+9 5E+00 

n ,W 5 q. ,a 

zcg g ^cg’ V 

0.2924E-01 

0.1223E-01 

0.1328E+09 

0.1200E+09 

0.1649E+02 

0.78U3E+01 

0.2331E+00 

WVW 

n swT 

0.7553E+00 

0 . 1+391E+00 

0.1259E+09 

0 . 3062E+09 

0.1727E+02 

0.1306E+02 

0.2938E+01 

n ,W ,q. ,a , 

zcg’ g u cg’ V* 
u 

1 1 

( 

0.1968E+00 

0 . 8606E-01 

0.1297E+09 

0.1576E+09 

0.1708E+02 

0.8976E+01 

0.7916E+00 

n ,W ,a , , 

zcg’ g’^cg’ v* 

0 

0.1+13 8E+02 

0.1730E+02 

0.1282E+09 

0.921+4E+10 

0.1233E+03 

0.281+6E+03 

0.1401E+03 

n ,W ,q , 

zcg g’ u cg’ u v’ 

BMRxlO -5 

0.1866E+02 

t 

0.8563E+01 

0.1199E+09 

0.1+008E+10 

.5860E+02 

.1180E+03 

. 5999E+02 


Cr\ 

00 














































































Table 6b Full State Feedback (Cruise Condition) 


Q 

R 

n zp (g>2 

COVARIANCE 

n zcg (e)2 

COVARIMCE 

2 

BMR(n-m) 

COVARIMCE 

S R (n-in) 2 

COVARIANCE 

W - 5 ^ 2 

COVARIANCE 

— 

<$ z (dg) 2 

COVARIANCE 

6 e (ag) 2 

COVARIANCE 

( 

01 

/.04 ^ 

.06 
• 25j 


.1821E-06 

•3781E-06 

•9973E+10 

.1297E+09 

.1483E+02 

.1381E+02 

. 3727E+00 


tk ) 

4T- 

CJ\ 

ro 

\ji 


.2972E-04 

.1588E-04 

.7078E+10 

.8506E+08 

.9280E+00 

.1778E+01 

. 2201E-01 
1 




SENSOR NOISE STUDY 
CRUISE CONDITION J = JY 


y = (n ,q ,a ,W ) T 
zcg eg’ v’ g 


Table 7a 


SENSOR NOISE STANDARD DEVIATIONS 


a q 


eg 


t-o/s 


a a : 
v 


.2 


o 


DESIGN 

TYPE 

J 2 

K SUCH THAT c 

? = -Ky 

6 e 

EIGENVALUES 

NO 

SENSOR 

NOISE 

IN 

OPTIMIZATION 

PROCEDURE 

168 

•966 -.163 1.31 -.266 

-4.09 2.28 -.854 2.68 

2.32 -1.53 .0656 -.00130 

~ 2 . 64 ±j 63.7 - 40.2 

-9.37*349.6 -2.26 

-1. 10+ j 36.0 -.0156 

- 6 . 83 * 331.2 - 13.0 

- 4 . 24 * 324.3 -.0000631 

- 21 . 7 * 314.0 

SENSOR 

NOISE 

ACCOUNTED 

FOR 

IN THE 

OPTIMIZATION 

PROCEDURE 

2 1+1 

3.07 0 . 0000461 - 0.0000462 0.311 

-1.67 -0.0000168 0.000420 2.25 

2.26 -0.000273 0.000236 -0.00904 

-1.71*365.2 -0. 00644+ j 0.0292 

-10.7*350.5 - 41.2 

- 2 . 85 * 338.4 - 26.5 

- 3 . 92 * 331.6 - 22.9 

-7-01±j21.0 
- 1 . 41 + 33 . 61 


-4 

H 




Table 7b 


^*^OVARIAN CE 
Design^^^ 

n zp p 
(g 2 ) 

n 

zcg p 
(g 2 ) 

BMR 0 

(n-m) 2 

V-) 2 

SsA (as) 2 

Sag) 2 

6 

e (dg) 2 

DESIGNED 

WITHOUT 

NOISE 

TESTED 

WITHOUT 

NOISE 

.1223E-03 

.2131E-03 

.8406E+10 

.1150E+09 

.1008E+00 

. 1625 E +01 

.3919E-02 

DESIGNED 

WITHOUT 

NOISE 

TESTED 

WITH 

NOISE 

.3226E+01 

. 6966 E+OO 

.112735+13 

.174735+11 

.6492E+02 

. 4869E+02 

.1414E+02 

DESIGNED 

WITH 

NOISE 

TESTED 

WITHOUT 

NOISE 

.2069E-03 

.2752E-03 

.7439E+10 

.9551E+08 

.343138-01 

. 5674E-01 

! 

.1499E+01 

DESIGNED 

WITH 

NOISE 

TESTED 

WITH 

NOISE 

.2070E-03 

.2752E-03 

.7439E+10 

.9551E+08 

.5674E-01 

. 1499 E +01 

.2515E-02 




- . 4321 1.087 16.00 

-6.484 -5.1*53 -29.il* 

5.500 -2.703 -45.73 

-3.47±j63.6 

-17.6±J43.6 

-2.00+J37.8 

-.000337±j26.4 

-5 .86±jl6 .6 

-. 00749± j .0933 

-40.9 

-31.3 

-10.5 

-.2251* 1.078 15.70 

-7.006 -5.1*91 -28.39 

5.1*61 -2.724 -45.38 

-3.07±j63.8 

-17 . 6± j 43 • 4 
-3.74±j46.5 

-.00210=026.4 

-5.97±ol8.5 

-.00752±o .0929 

-4l.O 

-31.3 

-10.5 






Table 
/ in 6 \ 

9a SENSITIVITY STUDY 


.966 

-.163 

1.31 -.266 

f -m 6 

m 

; y = (n , q , a . W ) 

» zcg’ u cg 5 v* g 

K = 

-4.09 

2.28 

-.854 2.68 

| 10 

\ 


2.32 

-1.53 

.0656 -.0013 


CRUISE GAINS EVALUATED AT VARIOUS FLIGHT CONDITIONS 


CONDITION 

cov 

n 

ZP 

COV 

n 

zcg 

EIGENVALUES 

CLIMB 

.3733E-03 

.3183E-03 

-2.58±j62.0 -.0122±j.0271 

-5.4o±j46.8 -4 o.o 

-2 . 8l± j 39 • 6 -6.68 

-2 .08-j32 .0 -1.08 

-3.l8±j21.7 

-23.7^3.36 

' 1 

CRUISE 

.1223E-03 

.2131E-03 

-2.64+J63.7 -40.2 

-9 • 37± j 49 . 6 -2.26 

-1.10±j 36.0 -.0156 

-6.83±j31.2 -13.0 

-4.24±j24.3 -.0000631 

-21.7±,3l4.0 

LAND 

.2012E-02 

.1145E-02 

-2.64+j6l.7 -4o.o 
-5.l4±j46.4 -24.4 
-2.75±j4o.3 -24.1 
-1.89+332.0 -4.30 
-3.09±j2l.6 -.983 
-.0127±j .0635 


[ 




Table 9b SENSITIVITY STUDY 


Q = 


10 


10 " 


K = 


-1.24 

2.94 


E(WWT) 



-.834 

-8.33 

3.50 

.706 

1,48 

1.94 

-.804 

-6.97 

2.62 


CRUISE GAINS WITH PROCESS NOISE 


FLIGHT 

CONDITION 

COVARIANCE 

n 

zp 

COVARIANCE 

n 

zeg 

X 

EIGENVALUES 

CLIMB 

0.5787E-02 

0.5748E-02 

-3.42+j6l.O -.000723±j .151 

-5.n±j46.6 -39.6 

-3.50±j40.8 - 28.7 

-3.56±j33.9 -21.1 

-2.48+J26.4 
-, 8 l 8 ±j 5 .58 

CRUISE 

0.2124E-03 

0.3388E-03 

-4.94±j6l.5 -.00747±j.0704 

-9.40±,j50.0 - 38.7 

-2,37±j38.4 

-8.80±j35.3 

-3.87±j22.7 

-23.1±J5.l6 

-1.68±J10.5 

LAND 

0.9490E-02 

0.8529E-02 

-3 . 26 ± j 60.9 - . 00117 ± 3 . 179 

-4.85±j46.0 -39-7 

~3.37±j4l.4 - 29.6 

-3.27±j33.7 -20.3 

-2 .37-^21 .4 
-.763±j4.63 




